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ABSTRACT 


On  August  20,  1978,  the  Synthetic  Aperture  Radar  (SAR) 
or.  board  the  satellite  SEAS  AT ,  gave  evidence  of  high  energy 
internal  oceanic  activity  at  the  shelfbreak  in  the  northern 
part  of  the  Bay  of  Biscay. 

Quantitative  spatial  measurements  of  internal  wave 
patterns  were  correlated  with  conventional,  quasi-synoptic 
in  situ  da-a  and  yielded  phase  speeds  of  0.551/S  for  high 
frequency,  nonlinear  internal  waves,  generated  in  groups  at 
the  canyons  indenting  the  shelfbreak  and  at  the  time  of  low 
tide  in  Brest.  Their  periods  were  found  to  be  of  70MIN  and 
their  amplitudes  estimated  to  be  15M.  Ocean  swell  refrac¬ 
tion,  observed  on  the  SAR  image,  together  with  a  localized 
significant  increase  in  wave  height  detected  by  -.he  SEASA? 
altimeter,  was  related  to  the  shear  of  a  northwestward  geos- 
trophic  jet  (ca.  3.5M/S)  at  ‘he  break. 

Possible  operational  consequences  of  these  features  are 
summarized  as  an  assessment  of  a  3AR's  capability  *0  give 
tactical  as  well  as  scientific  real-time  information  on  the 
internal  ocean. 
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RESUME 

Le  20  Acut  1978,  Is  SAR  (Synthetic  Aperture  Radar) 
embargue  sur  le  satellite  SEASAT,  a  sis  an  evidence  une 
intense  activite  interne  de  1* ocean  au  voisinage  du  talus 
continental  dans  la  partis  nord  du  Solfe  de  Gascogne. 

Des  assures  guantitati7 es  des  manifestations  de  surface 
des  ondes  internes  ont  ete  correlaes  avec  des  donnees 
conventionnelles  "in-situ”,  et  quasi  simultan^es.  Des 
vitesses  de  phase  de  0,55  metres  par  seconde  ont  ete  mises 
en  evidence  pour  ies  trains  d*  ondes  internes  non-linsaires, 
gene res  a  l'heure  de  la  basse  aer  de  Brest,  dans  les  canyons 
bordant  le  plateau  continental;  la  periode  de  ces  ondes  est 
de  70  minutes  et  leur  amplitude  a  ete  estia&e  a  15  metres. 
La  refraction  de  la  houle  observes  sur  1'  image  SAR,  en 


conjonction  avec  1*  observation  par  1‘aitimetre  de  SEASA? 
d' une  amplification  localises  mais  significative  de  la 
hauteur  de  vague,  a  ete  attribute  au  courant  de  cisaillement 
dun  "jet”  geostrophigue  de  surface  (da  l’ordre  de  0,5  metres 
par  seconde)  portant  vers  le  nord-ouest  le  long  du  talus. 
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&.  FOREWARD 


Remote  seasing  gives  satellites  an  essential  role  in  the 
monitoring  of  ocean  variability.  Synthetic  Aperture  Radar 
(SAR)  in  particular,  has  the  ability  to  display  images  of 
surface  manifestations  of  internal  processes  which  could 
typically  only  be  inferred  from  lengthy  and  costly  in-situ 
sampling  from  ships  at  sea  or  instrumented  buoys.  Sven 
though  the  dynamical  nonlinear  linkage  involved  in  ‘he  SAR 
imaging  mechanism  is  presently  incompletely  understood,  the 
vivid  and  spectacular  images  ;  e. g.r  as  those  acquired  by 
SEAS  AT,  give  a  tremendous  amount  of  information  about  well- 
known  or  totally  new  oceanic  phenomena.  To  study  such  an 
image  of  the  Bay  of  Biscay  (Pig.  1.1  and  1.2)  was  a  chal¬ 
lenging  and  exciting  topic  of  thesis  research  for  the  author 
who  has  experienced,  over  several  years,  the  sea-gcing 
approach  to  military  oceanography  in  this  area.  The  excep¬ 
tional  opportunity  for  doing  such  research  was  offered  at 
the  Naval  Postgraduate  School  where  direct  support  was  found 
and  the  introduction  to  other  facilities  was  given,  ^spe- 
cially  to  the  Fleet  Numerical  Oceanography  Center  in 
Monterey  and  to  the  Jet  Propulsion  Laboratory  in  Pasadena. 

B.  DATA  ACQUISITION 

The  optically  processed  SAR  image  was  the  starting  point 
of  the  total  data  set  collected.  It  was  forwarded  by  the 
Jet  Propulsion  Laboratory  (JPL)  of  the  California  Institute 
of  Technology  in  Pasadena,  and  a  first  heuristic  study  of 
the  image  gave  the  necessary  evidence  of  its  research  possi¬ 
bilities.  The  swath  is  100KN  wile  and  goes  rouahly  from 


Brest  in  Britan.ny  to  La  Coruna  in  Salicia;  internal  wav® 
activity  is  clearly  visible,  particularly  over  the  conti¬ 
nental  slope  (Pig.  1.2). 

Conventional,  quasi-synoptic  data  were  then  sought  and 
requested,  either  directly  through  the  French  Navy  authori¬ 
ties,  or  by  th®  Naval  Postgraduate  School  from  D.S.  and 
foreign,  military  and  civilian  institutions. 

A  fairly  dense  set  of  XBT's  (bathythermographs)  was 
collected  from  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  and  the  French  "Service  Fiydr  cgraphigue  ®t 
Oceanographique  de  la  Marine"  (SHOM)  .  The  format  of  these 
8T*s  gave  the  depth-temperature  pairs  at  turning  points; 
only  five  profiles  extend  below  533M.  Sixty  nine  XBT 
profiles  were  finally  considered,  corresponding  to  a  time 
period  of  seven  days  around  the  SAB  date  and  Located  in  the 
neighbourhood  of  the  image.  The  densest  spacing  covers  the 
continental  slope,  but  not  where  the  most  intense  activity 
is  noted  from  the  5  A  R  image. 

Local  atmospheric  data  were  obtained  from  the  French 
National  Weather  Service.  Global  weather  charts  over  the 
Atlantic  Ocean  and  swell  hindcast  for  the  Bay  of  Biscay  were 
obtained  from  the  FNOC  archive  center. 

Contemporaneous  moored  current  meter  time  series  were 
considered  as  important  auxiliary  information  regarding  the 
internal  baroclinic  modes  of  vertical  movements,  and  the 
mean  flow.  The  closest  mooring  was  a  subsurface  array  at 
1524,  3032,  3  976  M  depth  anchored  at  46<>57*2N,  009<>57'iw 

during  the  Northeastern  Atlantic  Deep  Sea  (NEADS)  mission. 
The  digital  da*a  (current  speed,  direction,  temperature, 
recorded  depth)  were  forwarded  oy  the  "Centre  National 
d' Exploitation  des  Oceans"  in  Brest,  where  the  "  Bureau 
National  des  Donnees  Oceaniques"  is  located. 


Other  satellite  sensors  cam  ad  by  SEASAT  acquired 
substantial  information  on  the  same  area  at  rapid  intervals. 
It  was  another  privilege  to  get  direct  computer  access  to 
the  Pilot  Ocean  Data  System  (PODS)  of  JPL,  thus  enabling  a 
thorough  survey  of  related  data  given  by  the  3  E AS  AT  a  1* im- 
eter  and  scatt  arose  ter. 

Climatological  data  were  blended  from  several  sources  in 
the  open  literature,  recent  papers  and  official  unclassified 
documents  of  the  O.S.  and  French  Navies. 

C.  OBJECTIVES 

The  main  aim  in  this  study  is  to  establish  a  definite 
correlation  between  the  SAR  surface  patterns  and  internal 
quasi-synoptic  features  established  by  other  means:  conven¬ 
tional  hydrographic  observations,  current  observations  and 
other  satellite  sensors. 

Avoiding  the  complex  problem  of  physical  linkage  in  the 
SAR  imaging  process,  a  pragmatic,  heuristic  approach  to  SAR 
analysis  was  made,  thus  tentatively  establishing  the  power 
of  SAR  to  monitor  the  oceans.  The  study  has  been  restricted 
to  the  upper  layer  of  the  Bay  of  Biscay  on  each  side  of  the 
continental  slope  and  shelfbreak,  where  both  the  SAR  image 
and  other  data  bases  were  found  to  provide  rich  information. 
Being  interested  only  in  internal  oceanic  features,  the 
swell  was  only  considered  in  regards  to  interaction  with 
quasi-steady  currents. 

In  the  course  of  the  research,  many  different  topics  had 
to  be  considered,  thus  greatly  enhancing  the  personal 
benefit  that  could  be  gathered  strictly  from  the  SAR  study. 
Conversely,  this  fact,  added  to  the  limirs  in  time  and 
scientific  expertise  of  the  author,  prevents  this  study  from 
b®ing  in  any  way  exhaustive.  Furthermore,  an  inherent 
feature  of  this  kind  of  quasi-synoptic  study  is  its 


study  is 


inability  to  support  objective  "sea-truth".  All  the  results 
presented  h-=re  must  be  thought  of  as  describing  a  par*icular 
event  at  a  particular  time  and  in  a  particular  environment 
that  can  only  be  considered  as  part  of  the  ever-growing 
information  available  to  oceanograpaers. 

Overall,  the  objective  of  a  thesis  is  no  increase  the 
knowledge  of  the  author  himself  and,  in  this  case,  in  topics 
related  both  to  the  use  of  conventional  methods  and  of  prom¬ 
ising  satellite-related  techniques.  Consequently,  this 
research  was  a  learning  experience  and  ar.  enriching  intro¬ 
duction  to  the  new  era  of  operational  military  satellite 
oceanography. 

D.  OBITS  ROTATIONS  CONTENT  IONS 

Onless  explicitly  stated  otherwise,  the  following  units 
will  be  used,  with  the  S.I.  system  as  a  general  basis. 

Quits  of  time 

GMT  Greenwich  time  (Z)  in  Hours,  Minutes,  Seconds. 

(The  local  summer  time  in  Prance  is:  GMT-020000; 
the  solar  time  at  the  SAR  image  location  is:  GMT-002000) 
Day-Hcur  groups  are  expressed  as  DDHHM1,  implicitly  days  of 
August  1978. 

Onits  of  length 
Depths  in  meters  (M) 

Distances  in  kilometers  (KM) 

Onits  of  t em oerature . 

Degrees  Celsius  (°C)  . 

Geographic  coord ir.ates 

Geographic  units  are  expressed  in  decimal  degre=s  of  lati¬ 
tude  and  longitude  (lat,lon).  Longitudes  are  counted  from  0 
to  360  in  most  computations. 
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An  auxiliary  Cartesian  grid  was  adjusted  over  a  nap  proiec- 
tion  of  the  area  to  meet  the  needs  of  the  particular  SAR 
image  orientation  and  general  bathymetric  gradient  (Pig. 
1.1)  This  grid  is  in  units  of  kilometers.  The  x-axis  is 
oriented  along  the  shelfbreak  with  positive  values  towards 
the  southeast  (directed  seaward  in  the  Bay)  . 

The  y-axis  is  oriented  along  the  center  of  the  SAR 
swath  and  *ovards  the  northeast  (directed  onto  the  French 
continental  shelf)  . 

Hence  x  values  give  the  relative  positions  on  each  side  of 
the  SAR  image  :(|x|<  60  for  locations  on  the  ground  track), 
y  values  give  relative  positions  on  each  side  of  the  shelf- 
break. 

The  transformation  was  established  as  follows  : 

Two  geographic  points  clearly  visible  on  the  SAR  image  and  a 
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ve  offshore 

bat  hymetric 

sounding  were  chosen  t 

get  the  linear  transform 

pa  rameters. 

la* 

47.80 

x*000 . 03 

Ion 

355.67 

y=1 32 . 42 

Penmate* h  point. 

la* 

43.59 

X=000 . 03 

Ion 

352.74 

Y=- 388.92 

Los  Cairos  point. 

la* 

46.00 

X=1  16.7 

Ion 

356.00 

Y=000 . 03 

Offshore  reference. 

The  corresponding  relationships  are  : 

lat  *-0.306261  X  ♦  3.308075  Y  *  U6.73 
Ion  =  0.009205  X  ♦  0.005620  Y  «•  354.92 


X  *-51.315  lat  +  7  3.73  Ion  -  23772. 


y  =  34.050  lat  ♦  57.17  lor. 


24217. 
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Starting  with  a  survey  of  the  recant 
SAR  principles  and  its  capabilities 
features,  the  various  assessments  tad®  of 
anisms  of  interaction  between  an  internal 
short  wind-driven  capillary/ultragravity 
can  detect  are  reviewed. 
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A.  SYITHETIC  AHTBHNAS 

Synthesizing  a  radar  antenna  consists  of  using  a  moving 
antenna  of  limited  size  and  combining  all  information  gath¬ 
ered  on  a  given  point  of  the  earth  during  a  given  time 
interval  (Fig.  2.2).  The  distant®  of  travel  on  the  orbit 
gives  the  length  of  a  virtual  antenna  of  much  finer  resolu¬ 
tion.  The  resulting  resolution  is  independent  of  range  and 
speed,  and  it  is  theoretically  half  the  real  antenna's  aper¬ 
ture.  In  the  range  (cross -velocity )  direction,  the  length 
of  coherent  pulses  determine  (in  tie  usual  wav)  the  corre¬ 
sponding  accuracy. 

In  the  case  of  SEASAT,  the  data  were  transmitted  in  real 
time  by  the  satellite  to  one  of  the  special  processing  earth 
stations.  The  final  record  is  either  a  computer-compatible 
tape  format,  or  a  film  conversion.  The  photographic  output 
is  displayed  in  four  parallel  subswaths  to  accoun*  for  a 
progressive  decrease  in  the  radar  response  intensity  with 
range  (Fig.  1.2).  Some  particular  images  over  limited  areas 


have  been  digitally  processed  with  a  valuable  enhancement  of 


The  extraction  of  a  two-dimensional  surface  image  spec¬ 
trum  is  done  either  by  an  optical  Fast  Fourier  Transform 
[Gonzales  1979  and  Shuchman  1978  %  or  by  numerical  methods 
[algorithm  of  Hasselman  1980].  Considerations  of  ‘he  rela¬ 
tive  qualities  of  these  processing  techniques  car.  be  found 
in  Mitchel  [  1982  ]  and  Tomiyasu  [  1978]. 

Charact erist ics  of  the  SEASAT  SAR  are  given  here  (Table 
I  and  Pig.  2.  1)  . 


TABLE  I 

SB A  SAT  SAR  CHARACTERISTICS 


satellite  altitude 

-  - — , 

303  KM 

satellite  velocity 

7.2  KM/S 

nfar-polar  orbit 

108  degrees 

view  angle  from  nadir 

23  degrees 

swath  width 

103  KM 

wavelength 

.  235  M 

frequency  L-band 

1275  GHZ 

real  antennae 

11x2  M 

synthetic  antenna 

IS  KM  1 

antenna  polarization 

99  l  linear  horiz. 

contrast  ratio 

9  d3 

radar  oain 

35  d8 

processing  resolution 

25  M 

oul3a  repetition  rate 

5  03  microseconds 

bulse  length 

33. microseconds 

OCEAN  IMAGING  MECHANISM  S 


The  application  of  SAR  systems  to  oceanography  requires 
extension  of  the  basic  SAR  earth  imaging  theory  to  account 
for  scattering  from  a  time-varying  surface  tha*  satisfies 
nonlinear  hydrodynamic  constraints.  3n  the  other  hand,  the 
general  scattering  models  for  low  resolution  radars  do  not 
apply  to  the  spaceborne  high-resolution  SAR's.  Th= 
complexity  of  the  problem  results  from  the  widely  nonlinear 
electromagnetic  and  hydrodynamic  interactions.  But  this 
nonlinearity  is  at  the  same  time  the  main  reason  for  the 


ror 


remarkable  capabilities  of  SAR  for  detection  of  mesoscale 
oceanic  phenomena. 

1.  Doppler  Effects  fi2A  lh°  22Ji23  Ocean  Sarface 

The  fact  that  the  ocean  surface  is  not  rigid  has  two 
effects:  ( 1|  defocusing  of  the  SAR  image  and  (2)  apparent 
displacement  of  the  radar  echo,  as  can  be  seen  by  comparing 
echos  of  moving  ships  and  their  stationnary  wakes  (Cf  para 

V-B-  1 )  . 

Empirical  studies  have  shown  that  to  recover  the 
optimal  focus  of  SAR  imagery  over  ooean  surfaces,  the  best 
speed  correction  to  use  is  the  long  wave  phase  velocity,  and 
not  as  would  seem  intuitively,  the  long  waves  orbital 
velocity  on  which  the  capillary  short  waves  ride.  In  the 
case  of  internal  waves,  this  phase  velocity  is  giver,  by 


where  "h"  is  the  depth  of  the  thermocline  and  "fa."  are 

the  densities  each  side  of  the  thermocline.  This  displace¬ 
ment  effect  will  periodically  modulate  the  image  intensity 
backscattered  by  long  waves,  causing  in  some  cases  periodic 
dark  and  light  bands  associated  with  the  pattern  of  the  real 
long  wave  pattern.  In  the  case  of  a  "closing"  current 
(cross-track  component  toward  the  satellite)  the  surface 
image  will  be  artificially  translated  in  the  positive 
along-track  direction;  hence  if  tais  "closing"  current  is 
ir.creasina  along-track  (fig  2.3),  the  value  of  this  transla¬ 
tion  will  increase  accordingly  and  there  will  be  stretching 
of  the  image  and  a  lower  Be agg-scat t er  than  expected. 


2.  Scatter  ir.g  of  Radar  Signal  from  Ocean  Wares 

Based  on  a  review  [Valenzuela  1978]  of  elsc-rcmag- 
r.etic  wave  scattering  from  ocean  waves,  the  standard 
approach  considers  the  ocean  surface  as  statistically 
described,  bat  a  deterministic  formulation  might  be  more 
appropriate  when  the  resolution  ar^a  is  small  enough  to 
correspond  to  a  few  wavelengths  only.  The  former  descrip¬ 
tion  is  nevertheless  easier  and  most  of  the  ^ime  more  accu¬ 
rate.  The  statistical  approach  may  lead  either  to  a 
Bracrg-rssonance-li  ke  mechanism  or  a  "facet''  mechanism 
(tangent-plane  theory)  [Brown  1976  ]. 

3.  Bragg  Model  vers  us  Tar. gent  Plane  Model 

The  Bragg  mechanism  is  limited  to  rough  surfaces 
satisfying  the  Rayleigh  criterion  : 


2.  -U  W  (rms'j  Cas  Q 
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where  "H"  is  the  typical  surface  height  fluctuation  and 


iiQh 


is  the  incident  angle  cf  the  radar  wave  on  the  water 


surface.  The  total  backscatter  is  then  given  by 


cr° (e».  ^  |g. .  (9) |  ty(<x 
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where  is  the  two-dimensional  Fourier  spectrum  of 

ths  sea  surface  with  an  ocean  wavenumber  "2kcos^”  and  "g-  •  " 

') 

represents  '■he  radar  reflectivity  of  ths  air-water  inter¬ 


face.  The  alternative  approach  is  the  tangent  plane  or 
facet  model  applicable  to  gently  undulating  rough  surfaces 
and  small  (§<<  30°)  incidence  angles  (for  SEASAT  SAR,  6  = 
23<>)  . 

We  have  in  this  case 
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where  R(0)  is  the  Fresnel  coefficient  for  normal  incidence 
and  ',<s2>"  is  the  focal  slope  variance  for  ocean  wav^s  of 
long  period.  In  both  cases  *he  intensify  variations  of  the 
SAR  signal  will  be  correlated  to  tie  energy  density  "H*"  of 
the  ocean  waves.  (They  will  be  at  the  same  time  modified  by 
the  mean  incidence  angle,  resulting  in  a  decreasing  average 
intensify  along  a  cross  track  line.)  Therefore,  the  nonli¬ 
near  periodic  energy  tra.nfer  from  the  capillary  waves  to  *:he 
mesoscale  phenomena  (currents,  swell  and  internal  waves) 
will  dramatically  affect  the  SAP.  imagery. 

C.  COUPLING  BITH  INTERNAL  OCEAN  DINAHICS 

Small  coherent  disturbances  in  the  current  field  can 
cause  important  modulations  in  the  ultragravity  wavefield 
[  Mollc-Christansen  1982]  :  e.g.,  for  wavelengths  of  30  cms, 

a  shear  current  of  0.003  S/sec  can  induce  a  modulation  at 
the  surface  to  which  a  SAP,  will  be  sensitive. 

1  •  'The  Sur  f  ace  Strain! no  atnd  the  3unchino  Effects 

The  possible  mechanism  involved  appear  to  be  either 
an  hydrodynamic  effect  (the  surface  straining)  ,  or  a  Doppler 
effect  (the  bunching)  . 

The  straining  effect  applies  to  all  infernal  ocean 
phenomena  (internal  waves,  shear  currents,  eddies). 
Variations  of  "  "  along  the  direction  of  wave  propagation 

are  induced  by  hydrodynamic  interactions  between  short  reso¬ 
nant  waves  ( A -  3W )  and  large  ocean  periodic  motions 
(X>100H).  This  mechanism  can  be  thought  of  as  an  alter¬ 
nating  pattern  of  divergence-convergence  of  the  free 
surface,  riding  on  a  long-wave  modulation  (sw=ll  or  internal 
wave)  .  The  resonant  waves  are  saturated  locally  at  the 
region  of  maximum  convergence  of  velocity,  ahead  of  the  wave 
crest  (Fig.  2.4)  . 


the  ito'iu la¬ 


in  the  short  wave  (Phillipsi  spectrum, 
tion  transfer  function  can  be  written  [Alpers, Ross 
Ruff  enach  1 981  ]: 


JlH 


a  n 


where  "|t**  is  a  time  constant  related  to  the  resonant  waves, 
is  the  angle  between  the  flight  (azimuth)  direction  and 
the  direction  of  propagation  of  the  internal  wave,  and  "Jl" 
is  the  angular  frequency  of  the  wave  component. 

The  velocity  bunching  affect  is  the  result  of  arti¬ 
ficial  displacement  of  the  image  of  periodically  moving 
surface  [Alpers  and  Rufenach  1979,  Alpers  et  al.  1991], 
Simply  put,  the  Doppler  shift  due  to  surface  velocity  is 
attributed  to  location  (azimuth)  and  the  sianal  is  mislo- 
cated  in  the  image.  It  seems  to  be  the  best  model  for 
azimuth  traveling  swell  and  eventually  for  internal  waves. 
However,  it  is  a  geometric  effect  that  does  not  account  for 
energy  transfer  as  well  as  *he  surface  straining  dynamic 
mechanism. 


2.  Applicat ion  to  Interna],  Waves 

The  general  form  of  for  a  wave  component  of 

wave  number  "K" ,  angular  frequency  "51",  and  amplitude  "A", 
may  therefore  be  written  more  generally  as  ' Vesecky  and 
Stewart  198  2]: 

<T8  =  <  <T*>[l+  A-  cos(K.r*-il-V^)]/ 

U  gcctJ  l*±]  y 

where  "3"  is  the  total  modulational  transfer  function  and 
can  be  considered  as  the  result  of  two  different  mechanisms 
(neglecting  the  tilting  effect  applicaole  only  to  swell): 

R  =  Rfhydrcdvn.)  ♦  R  (vel.  bunching) 
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of  short  capillary/gravity  waves  to  form  slicks,  and 
(2)  The  surface  current  field  induced  by  the  internal  waves 
result  in  a  periodic  modulation  of  the  height  of  the 
capillary/gravity  waves. 

ActuaLly  the  two  possible  mechanisms  implied  give 
contradictory  results:  damping  oy  oil  slicks  will  occur 

over  the  convergence  lines,  whereas  the  current  hydrodynamic 
damping  will  be  localized  on  the  divergence  areas.  It  is 
the  author's  assumption  that  whenever  the  surface  patterns 
visible  on  the  SAR  image  are  linear,  localized  and  filamen¬ 
tary  (Cf  patterns  described  in  7-B-1 i  ,  the  oil  slick  expla¬ 
nation  is  appropriate,  whereas  alternation  of  high  and  low 
brightness  can  be  more  easily  explained  by  hydrodyn amical 
damping  due  to  internal  waves  (Cf  group  pat-erns  described 
in  V-B-3  and  7-8-U)  . 

According  to  the  surface  straining  mechanism,  an 
internal  wave  traveling  at  phase  speed  "c"  produces  a 
surface  straining  and  a  velocity  pattern  "u(x-ct) "  moving 
with  the  waves.  Short  surface  gravity  waves  travelina  in 
the  same  general  direction  as  the  internal  waves  can 
strongly  interact,  particularly  when  the  group  velocity  of 
the  surface  wave  is  close  to  the  phase  velocity  of  th° 
internal  wave. 
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Gargett  and  Hughes  [1981]  give  another  descr ip-  icr 
the  same  principle.  The  horizontal  surface  currents 
associated  with  internal  waves  are  Large,  being  of  the  same 
order  as  the  wave  phase  velocities  themselves.  They  are 
oriented  approximately  either  parallel  (in  the  troughs) 
antiparallel  (on  the- crests),  to  the  direction  of  propaga¬ 
tion.  A  wind  field  superimposed  on  the  current  will  excite 
a  capillary/uL tragravity  wave  spectrum  that  is  steeper  where 
the  wind  opposes  the  internally  generated  current  than  where 
it  is  parallel  to  it  (Fig.  2.4).  Eventually,  a  modulational 
interaction  couples  the  surface-wave  spectrum  to  that  of  the 
internal  waves  and  possibly  creates  a  growth  of  the  latter 
[Dysthe  1981]. 

Indirect  detection  of  internal  waves  may  be  possible 
by  analyzing  the  variations  of  the  swell  energy  spectrum 
directly  extracted  by  FFT  of  the  SAR  data.  Refraction  of  the 
swell  is  then  the  mechanism  by  which  the  internal  waves  are 
manifesoed  on  the  sea  surface. 

D.  SE1SITIVITY  OF  &  SIR  PO  INTERNAL  MOTIONS 

1  •  Environments  1  Limits  t  i  on  s 

For  possible  radar  imaging  of  ocean  internal  wave 
motion,  the  first  essential  condition  is  the  existence  of 
surface  capillary  or  ultra -gravity  waves.  This  excludes 
total  calm  wind  conditions,  but  not  rough  and  even  breaking 
seas.  The  relative  directions  of  the  short  surface  waves 
and  th®  internal  wave  field  affects  the  hydrodynamical 
interaction  process.  Furthermore,  the  displacement  of  the 
satellite  may  induce  Doppler  imaging  effects  deper.deno  or. 
the  angle  which  the  radar  waves  make  with  the  surface  waves. 
Hence  ,  for  a  given  infernal  wave  field  and  known  satellite 
orbit  parameters,  the  optimal  wind  velocity  coalitions  could 
be  stated.  It  seems,  but  if  would  need  exDerimenfal  cor.fir- 
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mation,  that  a  wind  blowing  parallel  to  the  propacrat  ior.  of 
the  internal  waves  is  optimal  for  strong  interaction.  The 
radar  modulation  is  largest  for  internal  waves  traveling  in 
azimuthal  direction,  but  on  the  other  hand,  it  is  assumed 
that  a  satellite  track  parallel  to  the  same  internal  wave 
field  may  induce  enhancing  Doppler-ef feet s . 

That  these  geophysical  factors  are  sometimes  contra¬ 
dictory  does  not  seem,  however,  to  inhibit  the  capability  of 
a  SAR  to  display  surface  manifestations  (observed  in  many 
different  cases)  of  internal  motions.  In  particular  the 
analysis  conducted  in  chapter  V  shows  internal  waves  propa¬ 
gating  in  almost  all  relative  directions,  thus  proving  that 
at  least  some  interaction  and  radar  modulation  exists  for 
every  configuration. 

2.  Prooitious  Internal  Wav®  Pr o o er ties 


The  group  velocity  of  wind-driver,  short  waves  is 
by  i  a  , 

c3  =  Tk  +  krs 

where  "T"  is  the  surface  tension.  The  optimal  internal  wave 
phase  soeed  should  then  be  of  the  same  order  such  that 
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implying  approximately  for  ultra-gravity  waves 
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which,  in  the  case  of  the  30C1  wavelength  waves  that  STASAT 
SAR  can  detect  gives 

AP 


P  ^  4*H 


For  a  two-layer  ocean  with  a  thermocline  at  501  we  obtain 

A? 


n 


V.  O.S  \0*\ 


The  surface  horizontal  particle  velocities  must,  on  the 
other  hand,  be  large  enough.  3ut  we  know  they  are  directly 
related  to  the  amplitude  of  the  internal  wave  itself. 

He  therefore  think  that  detection  of  internal  waves 
by  their  surface  manifestation  should  be  favored  by: 
a  deep  mixed  layer  (H)  , 
a  sharp  theraocline,  and 
a  density  jump  of  the  order  of  W  43  a 
(H  in  meters) 

3.  Favorabla  Cases  :  an  Assessment 

Internal  tidally  generated  waves  on  the  continental 
shelf  have  frequency  properties  and  amplitudes  than  cause 
them  to  fall  info  the  stated  favorabla  conditions.  Those 
internal  waves  are  generated  most  commonly  in  the  highly 
stratified  seasonal  fhermocline  by  interaction  between  a 
barotropic  internal  tide  and  bathymetric  irregularities. 

As  a  particular  example,  the  continental  shelf 
offshore  France  in  the  Bay  of  Biscay  has  strong  tidal 
currents  relative  to  a  semidiurnal  tide  propagating  north¬ 
ward.  The  day  (20  August  1978)  of  passage  of  SSASAT  over 
*his  area  corresponds  to  a  spring  tide,  and  the  bathythermo¬ 
graphs  show  the  strong  stratification  common  to  the  summer 
season.  A  search  of  meteorological  data  has  confirmed 
average  summer  conditions  to  have  prevailed,  with  probable 
winds  from  the  northeast.  Previous  conventional  studies 
have  shown  chat  internal  waves  are  frequently  observed  and 
propagate  perpendicular  to  the  SH/SE  orientation  of  the 
continental  slope. 
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E.  COMBINED  UTILIZATION  OP  CONVENTIONAL  SENSOBS 


Adequate  knowledge  of  wind  conditions  is  necessary  for 
the  reliability  of  any  SAR  analysis.  Ship  observations,  as 
well  as  data  from  other  remote  sensors  (e.g.,  a  scat t erom- 
eter)  ,  were  used  in  this  study.  (Both  the  scat terometer  and 
SAR  are  side-looking  sensors.  Hence  they  offer  generally 
overlapping  coverage.)  The  required  knowledge  of  the  stra¬ 
tification  and  of  the  spectrum  of  periodic  movements  can  be 
derived  from  bathythermographs,  and  from  moored  instrumenta¬ 
tion  measuring  current  and  temperature  time-series. 

1 .  Bathythermographs 

Or.e  bathvthermographic  or  bathycelerimet ric  measure¬ 
ment  allows  the  computation  of  the  eigenvalues  solution  to 
the  internal  wave  generation  problem  '  Baines  1982  ].  It  is 
then  possible  to  find  the  wavelength  of  internal  waves  from 
the  spacing  of  wave-packets  generated  by  tidal  flow  over 
iopograohy. 

If  the  repetition  rate  of  the  forcing  "r"  is  known 
(tidal  internal  waves  for  instance)  one  can  derive  from  the 
spacing  between  packets  "<i"  ,  the  phase  speed 

Cp  ■=.  i  c 

A  bathythermograph  in  situ  makes  possible  the  computation  of 
the  main  density  gradient  on  the  thermocline  from 

Comparing  observed  and  theoretical  values  of  the  phase 
speed,  it  is  possible  to  infer  the  tnermocline  depth  by 


and  eventually,  by  combining  thi3  information  wi*h  sea- 
surface  temperatures,  the  heat  storage  in  the  ocean 
[ Mollc-Chri stens en  and  flascarenhas  1979  ]. 

2.  Current  and  Tempera  tare  Observations 

Time  series  of  current  and/or  temperature  measure¬ 
ments  at  several  depths  make  possible  the  vertical  modal 
analysis  of  the  vertical  and  horizontal  movements  at  a 
point , 

High  resolution  observations  of  the  sea  surface 
temperature  by  airborne  or  spacebone  infrared  radiometers 
might  detect  horizontal  gradients  due  to  variations  of  the 
surface  thermal  layer  depth  directly  related  to  the 
stretching  effect  of  an  internal  wave  field. 
Quasi-simultaneous  temperature  transacts,  in  the  sam»  way, 
can  resolve  the  space  characteristics  of  an  internal  wave 
field  that  can  only  be  inferred  from  a  fixed  observation 
station  together  with  internal  wave  theory. 


Pigure  2.3  D0PPL2B  2FP2CTS  3?  003825?  3321H. 


III.  GENERAL  GEOPHYSICAL  FEATURES  OP  THE  BAY  OP  BISCAY 

According  to  Sverdrup  et  al.  '1942]  ’* t h «=  terminal 

branches  cf  tie  Gulf  Stream  system  are  not  all  well  known; 
irregular  branches  Which  turn  to  the  south  terminate  in 
great  whirls  off  the  European  coast'*.  These  irregular  it  ies 
could  correspond  to  semi-permanent  circulation  patterns  or 
propagating  eddies.  Recent  studies  have  be  =  n  made  of  the 
mean  circulation  and  of  mesoscale  features  in  the  Bay  of 
Biscay,  particularly  in  the  vicinity  of  the  Armorican  conti¬ 
nental  slope,  from  which  some  general  conclusions  can  be 
drawn . 

Several  general  geophysical  features  relevant  to  the 
subarea  covered  by  the  SAP  image  will  be  consid-red,  giving 
particular  emphasis  to  the  summer  climatological  conditions 
under  which  the  SAR  image  was  taken. 


A.  ATMOSPHERIC  CIRCULATION 

Atmospheric  conditions  over  the  Bay  of  Biscay  are  deter¬ 
mined  by  the  location  of  the  North  Atlantic  Subtropical  High 
(pressure  center)  centered  in  February  around  25°N  and  in 
August  around  35°N.  The  High  extends  over  29  to  25  degrees 
of  latitude  in  summer  (vs.  only  15  in  winter)  and  is  limited 
to  the  north  by  the  Icelandic  Low.  The  general  atmospheric 
circulation  is  consequently  a  mean  eastward  flow,  weak  in 
summer,  which  is  modified  by  the  presence  of  the  European 
continent  and  perturbed  by  the  passage  of  fronts  originated 
ever  the  Nortwest  Atlantic  in  the  vicinity  of  the  Polar 
Front  (between  U5°and  50°N  in  summer).  Table  II  cives  the 
statistical  averages  for  August  an  t hr e=  stations,  respec¬ 
tively  north,  east,  and  southeast  of  the  3av  of  3iscav. 
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B.  BATHYHETRI 


The  Bay  cf  Biscay  (Gclfe  ds  Gascogne)  is  a  fairly  large 
bay.  If  is  450KM  (250NH)  wide  at  its  western  opening  on  the 
West  European  Basin  and  extends  553KM  (300  NM)  to  the  south¬ 
east  (Fig.  3.1).  Its  physiographic  characteristics  have 
been  presented  by  Laughton  [  1975  ]. 

The  abyssal  plain  has  an  average  depth  of  4500M, 
increasing  southward  to  49003.  It  is  limited  to  the  south 
by  the  narrow  (40KJ1)  Cantabrian  continental  shelf,  and  to 
the  north  by  a  wide  (200KB)  and  shallow  (130M)  shelf  (Fig. 
3.2) .  This  latter  shelf  extends  far  beyond  the  Bay  of 
Biscay  to  encoapass  the  British  Isles.  The  designation  of 
Celtic  Sea  is  generally  given  to  the  northern  part  of  the 
shelf;  the  Araorican  Continental  Shelf  is  to  the  south  of 
48® N,  the  latitude  of  Brest. 

The  "La  Chapelle  Bank",  at  tie  above  approximate  limit 
(48°N,  007ow>  is  100  to  1503  deep  and  characterized  by  large 
scale  and  smooth  ridges  (amplitude  ca.  40B)  oriented  yNB-SSW 
and  20K3  apart.  Its  southwestard  liiit  (4704l'N,  007013'  W)  , 
sometimes  referred  to  as  the  bank  per  se,  exhibits  sand 
waves  [Cartwright  1958],  correlated  with  scattering  layer 
undulations  [Stride  1960  ]  as  a  consequence  of  internal  tidal 
wave  activity  ’Carruthers  1  963  and  Cartwright  1963  ].  These 
sand-waves  have  an  average  height  of  9M  and  a  wavelength  of 
850B;  they  are  oriented  toward  113°r  (parallel  to  the  edge 
of  the  shelf)  and  extend  over  at  lease  20KS  from  the  slope. 
Southeast  of  La  Chapelle  Sank,  the  continental  seafloor  has 
a  variable  geological  nature  (rocks,  mud,  gravel  and  sand)  , 
homogeneous  in  depth  with  a  very  gentle  slope  (0.053S). 

The  shelfbreak  of  the  Armorican  Shelf  is  oriented  130°T 
and  occurs  between  180  and  2003  (Fig.  3.3)  The  wid^h  of  the 
continental  slope  decreases  from  35K3  at  La  Chapelle  Bank  to 
20K3  at  47<>N,  0Q5°30'W,  where  the  steepness  is  maximum  (2(H 


between  the  233  and  1000W  isobaths!  .  This  steep  slope  is 
indented  by  many  canyons  and  ridges  normal  to  its  axis.  The 
designations  given  by  Laughton  [1975]  will  be  used  here.  On 
the  swath  covered  by  the  SAR  image,  five  major  canyons  are 
visible:  the  Renmarc’h  Canyon,  tae  3uilvir.ec  Canyon,  the 
Odat  Canyon,  the  Blavet  Canyon  and  the  Belle-Ile  Canyon 
(Pig.  3.3). 

The  continental  rise,  south  of  the  slope,  is  200KM  wide. 
It  is  broken  by  the  Gascony  Knoll  at  45°20'N,  005°20'»,  with 
a  minimum  depth  of  3950M,  at  the  eastern  edge  of  the  SAR 
image.  The  Santabrian  Knoll  centered  at  45<>05'N,  008°00'W 

is  steeper  and  higher  (30Q3M  minimum  depth,  20  by  40KM  hori¬ 
zontally)  ,  but  it  is  50KS  off  the  western  edge  of  th9  SAR 
image. 

C.  OCEAN  CLIMATOLOGY 

The  surface  circulation  in  the  Northeastern  Atlantic  is 
not  as  well  defined  as  that  in  tae  western  (Gulf  Stream) 
regime  of  the  North  Atlantic  anticyclonic  circulation.  After 
crossing  the  mid-Atlantic  ridge,  tae  Gulf  Stream  loses  its 
coherency  and  splits  into  several  branches,  which  are  weak 
and  greatly  variable  on  all  space  and  time  scales  and  which 
eventually  decay  into  fields  of  poorly  defined  eddies  (Fig. 
3.4).  The  easternmost  of  these  branches  flows  south  along 
the  European  coast,  probably  creating  a  weak  cyclonic  circu¬ 
lation  in  the  Bay  of  Biscay,  as  it  weakens  along  the  north 
coast  of  Spain  [Millet  1980].  A  transient  westward  surface 
current  in  the  northern  part  of  the  3ay  foil  ws  long  periods 
of  eastward  winds;  if  is  referred  to  as  ...a  Rennel  Current 
and  described  in  the  French  "Instructions  "antiques"  [SHOM 
1975  ]  as  being  probably  narrow,  but  strong  (O.b  to  0.7^/S). 


The  ne ar-b3t*.oi  circulation  is  highly  variable. 
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However,  •‘■here  is  a  quasi  permanent  cyclonic  bottom  circula¬ 
tion  on  the  abyssal  plain  csnterel  at  44030'N,  004030  ''*, 

which  is  driven  by  the  eastward  flow  along  the  Iberian 
Peninsula.  Dn  the  Araorican  Shelf  the  mean  bottom  current 
drifts  to  the  northwest  from  March  to  July,  and  then  towards 
the  southeast  from  August  to  November.  The  average  current 
speed  measured  at  1M  above  the  sea  floor  is  0.04M/S.  [Camus 
1971  ]. 

The  intermediate  waters  in  the  abyssal  plain  of  the  Bay 
are  under  the  influence  of  the  Mediterranean  outflow,  which 
spreads  north  out  cf  the  Straits  of  Gibraltar.  This 
intruding  water  mass  is  cantered  at  a  depth  of  about  1100M 
and  its  T-S  characteristics ,  originally  T=11.9,  s  =  36.5Qo/oo, 
decrease  by  mixing  to  T=10,  S=35.33»/oo  in  the  greater  part 

of  the  Bay  of  Biscay.  Assuming  a  level  of  no-motion  at 
2009H,  results  from  two  separate  hydrographic  data  sets 
(July  1963  and  November  1957)  exhibit  the  following  charac¬ 
teristics  of  the  Mediterranean  Watar  circulation  [ Le  Floch 
1969  ]:  a  permanent  cyclonic  circulation  occurs  in  the 

southwest  portion  of  the  Bay  cf  Biscay,  approximately 
centered  at  44°30'N,  007O30,r4  on  which  is  superimposed  vari¬ 
able  meanders  and  probable  eddies;  the  speed  of  this  circu¬ 
lation  is  of  the  order  of  3.05M/S. 

D.  WATER  HASSES,  FRONTAL  STRUCTURES,  AND  EDDIES 

Hydrographically ,  the  deep  ocean  wafers  of  the  Bay  of 
Biscay  have  three  main  components;  the  relatively  warm 
surface  wafers  are  from  the  North  Atlantic  anficyclonic 
gyre;  the  deep  and  bottom  waters  are  derived  from  Antarctic 
Bottom  and  Intermediate  Wafer;  tne  midcoluan  water  mass 
charcterisf  ics  vary  with  the  degree  of  mixing  of 
Mediterranean  Wafers  with  either  one  of  the  first  two 
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components.  On  the  Armorican  Shale  a  mass  of  colder, 
saline  water  occurs  as  a  bottom  layer  adjacent  *.c 
French  coast  (Fig.  3.5). 

We  will  consider  as  a  front,  the  locus  of  sharp  hori¬ 
zontal  gradients  in  temperature  or  salinity  at  an  arbitrary 
depth.  Two  main  surface  fronts  are  driven  by  tidal 
currents,  one  offshore  Britanny  due  to  a  horizontal  decrease 
in  mixing,  and  the  other  at  the  continental  slope  which  will 
be  discussed  in  more  detail.  Two  pelagic  fronts  have  beer, 

identified  :  (1)  the  the  rmalhaline  front  delimiting  the 

extension  of  Mediterranean  Waters,  which  is  difficult  to 
monitor  but  which  may  have  important  acoustical  effects,  and 
(2)  the  thermal  front  bounding  the  coastal  bottom  cold  water 
mass  of  the  Armorican  Shelf,  which  is  enhanced  in  summer  by 
surface  heating,  but  which  has  no  surface  temperature  mani¬ 
festation  in  that  season. 

Eddies  occur  in  the  Bay  of  Biscay  in  summer  conditions, 
when  the  thermal  structure  of  the  upper  water  column  is 
strongly  developed.  Based  on  infrarad  satellite  images  of 
eddies  bordering  the  Ushant  frontal  zone,  these  eddies  are 
most  apparent  along  geostrophically  balanced  fronts  where 
cross  frontal  mixing  takes  the  form  of  cyclonic  disturbances 
[Pingree  1978].  They  have  a  length  scale  of  20  to  40  KM, 
with  a  time  scale  of  a  few  days  and  translational  velocity 
less  than  2CM/5.  In  the  case  of  the  Armorican  shelfbreax, 
irregular  eddies  form  on  each  side  of  the  230M  isobath;  in 
this  process  "the  cooler  water  sinks  below  the  warmer  water 
and  is  stretched  vertically'*  [Pingree  1978  ], 

Over  the  abyssal  plain,  less  apparent  eddies  have  beer, 
detected  by  drifting-buoy  experiments  and  by  satellite 
imagery  (Diskson  1981).  These  eddies,  both  anticyclor.ic  and 
cyclonic,  seem  to  partly  decouple  salinity  and  temperature 
features  from  density  surfaces.  From  an  AVHRR  image  taken 
on  14  May  1979,  and  in  the  v icini : /  of  the  SAP  image,  it 
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appears  that  cyclonic  eddies  were  aligned  or  continental 
rise  east  of  006°W,  and  that  a  pair  of  cyclonic  and  ar.ticy- 
cloric  eddies  was  located,  respectively,  northwest  and 
southeast  of  '45O30'N,  Q07O30»Sf.  The  length  scale  of  these 

deep-water  eddies  is  50  to  100KM  [Dickson  1981], 

E.  TIDES 

"The  tides  of  the  Atlantic  Ocean  can  be  considered  as 
composed  of  two  parts,  a  co-csciliat ing  tide  which  is  main¬ 
tained  by  the  tide  of  the  Antarctic  Ocean,  and  a  free  tide 
which  is  maintained  by  the  direct  effect  of  the  fide 
producing  forces.  As  a  result,  the  semidiurnal  fide  can  be 
considered  as  composed  of  two  standing  oscillations  having  a 
phase  difference  of  one  quarter  of  three  lunar  hours." 
[Sverdrup  1942]  The  amphidromic  point  is  "estimated  reli¬ 
ably"  [Cartwright  1980]  at  50°N,  339<>W,  so  that  the  fide 

will  propagate  from  south  to  north  in  the  Bay  of  Biscay  with 
a  90  degrees  Greenwich  phase  lag  in  the  south  versus  110  in 
the  north  (Fig.  3.6).  The  tidal  amplitudes  increase  from 
1.0M  offshore,  to  1 . 2t!  over  the  continental  rise  and  1.4M  at 
50KM  off  the  French  coast.  By  way  of  comparison,  the  diurnal 
tidal  amplitudes  over  the  area  are  of  the  order  of  0. 07* 
(Table  III) .  The  tidal  energy  dissipation  computed  by 
Cartwright  [  1980  ]  largly  exceeds  direct  estimates  based  on 
bottom  friction,  leading  him  to  conclude  that  "a  substantial 
loss  must  occur  along  the  Biscay  shelf  south  of  Ouessanf 
(Ushant  Island)."  And  indeed,  Pingree  [  1982]  infers  from 
both  observations  and  nonlinear  tidal  models  "-hat  "about  60^ 
of  the  tidal  energy  flux  onto  the  shelf  occurs  in  the 
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j  Greenwich  phase  lag  in  degrees. 
I  Tidal  amplitudes  in  CM. 


Internal  tidal  modes  have  been  studied  from  current 


meter  observations  at  47047N  010°13'ii  in  April  1978,  [James 
1982  ].  Results  support  the  notion  that  internal  tides  are 
dominated  by  the  lowest  vertical  modes  and  that  the  M2 
component  is  significantly  coherent  with  the  equilibrium 
fid®.  Barotro  pic  tidal  current  ellipses  were  found  to  be 
narrow  and  directed  perpendicular  to  the  shelf,  thus 
"forcing  the  isopycnal  surfaces  over  topographic  features 
(the  accepted  explanation  for  the  generation  of  internal 
tides)"  [James  1  992  ].  They  show  a  regular  turning  (counter¬ 
clockwise)  of  their  orientation  axis  from  570  to  3670M  depth 
and  could  thus  produce  baroclinic  activity.  The  tidal 
ellipses  are  themselves  of  clockwise  rotation  at  depths  less 


than  1270M. 


F.  IHTEHNAL  RIVE  MOTIONS 


Internal  waves  are  a  common  phenomenon  in  the  stratified 
ocean  and  exist  in  the  frequency  band  from  the  local  iner¬ 
tial  frequency  (16. 4H  period  at  483N>  to  the  maximum  Brur.t- 
Vaisala  frequency.  At  tijal  frequencies  (s.g.  .  such  as, 
semidiurnal),  periodic  motions  should  b*  called 


internal  (barocliaic)  tiles  as  they  result  from  forcing  by 
tidal  forces.  At  frequencies  greater  than  th*  Brunt -Vais ala 
frequency,  perturbations  cannot  propagate  as  waves  and  are 
turbulent. 

In  the  center  of  the  Bay  of  Biscay,  temperature  time 
series  f Hecht  and  Hughes,  1971]  have  shown  chat  semidiurnal 
internal  tides  are  greatly  enhanced  by  spring  (barotrcpic) 
tides  and  can  generate  high  frequency  oscillations  with 
marked  peaks  at  periods  of  T=7 1  and  2 4 :d I N .  These  deep-water 
internal  tides  could  be  generated  by  the  incidence  of  the 
barotrcpic  tide  on  the  continental  shelf,  to  then  propagate 
offshore  normal  to  the  coast  line  -'Rattray  i960  ].  Their 
amplitudes  would  then  be  given  by: 

4-^t- / 

where 

'  t  ii 

l  _  x\  -  \ 

and  "A"  is  the  amplitude  of  the  internal  wav  =  ,  " 2 H "  is  the 

tidal  range,  '*h"  is  the  depth  to  ‘■.he  thermocline ,  "D"  is  the 
deep-water  depth,  "dM  is  the  depth  of  the  shelf  (assumed 
constant),  "W*  is  the  width  of  the  continental  shelf,  "r."  is 
the  frequency  of  the  internal  wave,  "f"  is  the  Coriolis 
parameter,  ”3"  is  the  acceleration  due  to  gravity,  and 
"  ^  "  is  the  relative  density  difference  between  th°  two 
layers.  Hecht's  [1971]  observed  ampli-udes  exceed  ‘rhcse 
estimated  by  the  above  formula  (23H  obseved  vs.  15.8M 
computed  at  spring  tide).  On  the  continental  shelf  (U6o»j 
003°50W),  r.earbottom  current  records  occasionally  exhibit 
bursts  of  periodic  oscillations  (ca.  3 0 M I N )  followed  durina 
a  few  hours  by  apparent  whi*-.e  noise  fluctuations  [  Le  Floch 
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IV.  GEQPHISIC&L  SITUATION  DN  AOGOST  20  1978 

A.  HETE0H0L03  Y  AND  NIND  ST  BESS 

On  the  23th  of  August  1978,  a  high  pressure  ridge 
extended  from  the  Azores  to  Western  Europe.  The  average 
surface  pressure  over  the  Say  of  Biscay  was  1023MB  (Fig. 
4.  1)  .  A  strong  storm  generated  on  the  15th  of  August 
offshore  Labrador  was  dying  over  Iceland  (Fig.  4.3).  A  weak 
cold  front  was  300NM  west  of  the  Bay  of  Biscay  on  the  19th 
at  1200  GMT  and  probably  underwent  frontolysis  in  the  area 
of  interest. 

From  ship  reports,  the  evolution  of  the  in-situ  atmos¬ 
pheric  situation  (ca.  470N,005°W)  can  be  analysed  (Table 
IV).  Between  200300Z  and  200600Z  the  cold  front  was  present 
near  the  location  of  observations  at  46°4GN,  304°3QW. 

TABLE  IV 

METEOROLOGICAL  OBSERVATIONS  ca.47H,005W 


I 

I 

From  synoptic  charts  and  these  ship  reports  a  comparison 
I  can  be  made  to  select  the  SEASAT  sea tteromat ?r  winds  a*  the 

f  time  of  the  SA 8  image  (fig.  4.2).  At  202142Z,  a  local  high 

pressure  center  was  in  the  vicinity  of  480N,  010o  w,  '■has 

establishing  an  anticyclonic  surface  flow  decreasing  from 
!  c a.  10M/S  south  of  Penmarc'h  (north  of  the  SAR  image)  to 

!  ca.  5  M/S  at  45°N  006°W  (south  of  tie  SAR  image)  .  The  best 

estimate  for  the  direction  is  a  flow  towards  20QOT  (winds 
norr h-nor*heast) . 

^  B.  SHELL 

A  heavy  (mean  significant  wave  height  of  2.5M),  and  long 
(mean  wavelength  of  230M)  swell  had  originated  2000KM  away 
l  (Fig*  4.4)  and  was  currently  prevailing  over  the  sea  state. 

It  was  propagating  roughly  130°T  (from  shipboard  meteorolo¬ 
gical  reports)  .  The  generation  areas  were  estimated  by 
atmospheric  charts  (Fig.  4.3)  to  be  in  the  southwest  quad¬ 
rant  of  the  surface  pressure  low.  Distances  and  time  of 
propagation  with  respect  to  the  SAR  image  location  and  time 
yielded  estimates  of  the  wavelengths  that  could  be  expected 
in  the  Bay  of  Biscay,  August  20,  as  a  function  of  mean 
direction.  The  deep  water  surface  gravity  waves  equations 
used  are  the  following:  _ 

the  period  is  expressed  as  ~T~ -  \ j  jLH_ 

3 

JX 

2.TT 

thus  yielding  \  10 2 

9  T:* 

with  "DH ,  distance  from  the  generation  area  ar.d  time 

elapsed  since  the  original  wind  generation. 

Mo  swell  spectra  were  available  to  compar®  with  the 
patterns  visible  on  ■‘■he  SAR  image,  but  a  correlation  with 
significant  wave  heights  as  measured  by  the  altimeter  of 


the  group  speed  is  given  by 


--fv/ 


■3  z 


SEASAT  was  ooesible, 


The  most  interesting 


in  the 


altimeter  pass  of  210721Z  (ascending  node  passing  by  U6°N, 
005°W)  was  a  sadden  increase  in  significant  wave  heigh*  from 
an  average  of  2a  to  more  than  *1  over  a  short  distance 
(15KB)  at  position  47oi5«ji,  006°00*ii  (Fig.  4.5),  i.e.  near 
the  shelfbreak.  This  astonishingly  abrupt  variation  could 
hardly  be  related  to  any  atmospheric  mesoscale  process.  If 
will  be  shown  (Cf  71-rE)  that  it  may  be  the  consequence  of  a 
strong  interaction  with  current  shear  occurring  at  the 
shelfbreak.  According  to  sailers  experienced  with  this  area 
[Chiari  1982],  crossing  the  200M  isobath  is  often  accompa¬ 
nied  by  a  characteristic  change  in  sea-state.  An  historical 
reference  is  given  by  Le  Danois  [1943],  who  attributes  the 
denomination  of  La  Chapelle  Bank  to  an  observation  mad®  by 
the  "Sieur  de  La  Chapelle'*  in  1695  of  surface  turbulence 
wrongly  associated  with  extreme  shoaling.  This  apparent 
shoal  was  resighted  in  176  4  and  a  rock  was  marked  on  the 
charts  for  the  following  decades.  A  well  known  example  of 
such  a  current-swell  interaction  is  the  occurrence  of 
"giant”  waves  [Huthnance  1980  ]  near  the  shelf  edge  off  the 
East  coast  of  South  Africa.  These  are  thought  to  arise  at  a 
caustic  associated  with  the  refraction  of  surface  waves  by 
the  Agulhas  Current.  Ir.  the  same  way,  an  alongshore  current 
at  the  shelfbreak  in  the  Celtic  shelf  could  explain  rhe 
increase  in  wave  height  observed  there. 


C.  HEAR  ARD  TIDAL  CURRENTS 

Data  available  for  direct  estimation  of  the  mean  curren* 
and  its  variations  over  the  area,  are  provided  by  the  ieeD 
water  subsurface  mooring  at  46°57*!f,  009°57*  iw  with  current 

meters  at  1524,  3032  and  3976B  (350KM  off  rhe  SAP  ground 

track  and  180KB  offshore  of  the  shelfbreak)  (Cf  1-B).  Low 
frequency  currents  at  the  three  levels  do  not  exhibi*  any 
obvious  correlation  and  vary  over  ail  directions  wi-h  a 


time-scale  of  thrae  months.  The  naan  currants  in  the  Last 
10  days  before  tha  20th  of  August  ara  r aspect ively : 

41  15241  toward  010°T  at  0.051/S 
at  30  321  toward  270°r  at  O.OOJM/S 

(veering  330°T  in  last  3  days) 

It  3976(!  toward  315°r  at  0.0011/S 
On  a  shorter  time  scale,  progressiva  vector  diagrams  show  a 
semidiurnal  pariodicity  at  each  of  the  three  levels  (Fig. 
4.6)  ;  the  mean  speeds  are  : 
at.  15241  0.081/S 
at  3032M  0 .05 M/S 
at  3976M  0.0015M/S 

At  20  21  40,  time  of  the  SEASAT  passage,  the  instantaneous 
current  vectors  were  : 

at  15241  toward  315°T  at  0.081/S 
at  30  321  toward  185<>P  at  0.061/S 
41  39761  toward  270°T  at  0.0015M/S 
The  passage  of  SEASAT  on  20  August  followed  a  spring  tide 
ranging  110  on  the  Fra  no  h  coefficient  scale  [SIOM  1968  ] 
(i.e.,  amplitudes  1 55S  greater  than  the  average  spring  tides 
of  coefficient  95)  .  At  the  time  of  the  imag*,  the  surface 
tide  at  the  shelfbreak  was  at  high  tide  plus  four  and  a  half 
hours  (Table  V).  The  surface  currents  were  estimated 
towards  230  °T  at  0.85M/S  (Table  4.5),  [SHOM  1968], 

D.  SPATIAL  TBEH1AL  STRUCTURE 

The  thermal  structure  over  the  area  was  established  from 
a  set  of  69  XBT *  5  (Cf  1-B)  .  The  time  span  of  this  tempera¬ 
ture  sampling  was  from  0900  17  August  to  0500  24  Auoust. 

The  spatial  Limits  are  between  43°  and  u3°Norih,  and  10° 
West  and  th*  coastline.  The  distribution  of  this  sampling 
is  conveniently  cantered  around  the  radar  image,  in  partic¬ 
ular  in  th«  vicir.itv  of  the  continental  slcoe  and  with  a 
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maximum  density  on  the  eastern  half  of  the  ground  coverage, 
as  shown  by  a  plot  of  the  spatial-temporal  distribution 
(Fig.  4.7)  with  respect  to  the  distance  from  the  shelfbrealc 
(y-coordinate  defined  in  1-C  and  Fig.  1.1). 

1  •  Analysis  of  the  dean 

A  smoothed  transect  was  established  by  averaging 
along  the  x-axis  with  a  wy"  step  increment  of  20KM  (Fig. 
4.8)  .  The  isotherms  generally  sloped  down  from  the  shelf 
region  to  the  open  sea  (on  the  average  10M  over  100Kd).  A 
marked  mixing  was  apparent  in  the  region  of  th -  shelfbreak 
(-10<Y<10)  with  shoaling  of  the  isothermal  layers  above  the 
thermocline  (14°C)  and  sinking  below  it.  Dynamic  heights 
relative  to  the  smoothed  tr  msec4:,  yield  a  mean  surface  gecs- 
trophic  flow  relative  to  150d  of  the  order  of  .2M/S  and 
about  2QKM  wide,  towards  the  northwest  at  the  shelfbreak. 


TABLE  TI 

SURFACE  TIDAL  COBRENTS  AT  H7°24«N  006O05'« 


—  .  ^ 

1  Hour 

1 

Direction 

' 

Speeds  «/S 

•1 

_J 

i  GMT 

o  T 

Spring  i 

Naan  j  AUGUST  20 

L  ‘ 

1 

i  04 

45.  . . 

— 

“1 

r 

.  H igh  Tide 

05 

05 

120 

0.  25 

0.15 

|  0.  3 

06 

05 

190 

0.  4 

0.2 

1  0.45 

07 

05 

20  0 

0.  5 

0.  25 

|  0.55 

OB 

05 

210 

0.  6 

0.3 

|  0.65 

09 

05 

22  5 

0.  75 

0.4 

1  0.85 

10 

05 

23  8 

0.  75 

0 . 4 

|  0.85 

1  1 

05 

236 

0.  5 

0.25 

|  0.55 

1  1 

26 

25  0 

0.  4 

0. 2 

1  0.45 

1  2 

26 

29  5 

0.  2 

0  .  1 

1  0.2 

1  3 

26 

10 

0.  3 

0.  15 

|  0.3 

14 

26 

35 

0.  5 

0.25 

1  0.55 

15 

26 

50 

0.  5 

0.25 

1  0.55 

16 

26 

70 

0.  45 

0.25 

I  3.5 

1  7 

()  f\  T  * 

. . . H iah  Tide 

17 

26 

12  0 

0.  25 

0.15 

1  0.'3 

18 

26 

190 

0.  4 

3.2 

1  0.45 

19 

26 

20  0 

0.  5 

0.25 

I  0.55 

20 

26 

210 

0.  6 

0.3 

I  0.65 

2  1 

26 

22  5 

0.  75 

0.4 

I  0. 85 

I 

SAR 

image 

--  at  214i 

—  ■‘■o  wards 

2300  o.85 

M/S - 

1 

22 

26 

238 

0.  75 

0.4 

|  0.85 

23 

26 

23  6 

0.  5 

0.25 

1  0.55 

23 

47 

250 

0.  4 

0.2 

1 _ 0-u 

_j 

The  geometric  approximation  of  the  shelf  by  a 
straight  line  in  the  averaging  proress,  and  probable  varia¬ 
tions  in  time  and  space  of  tnis  mixing  region,  reduce  the 
more  dramatic  features  that  _^.n  be  seen  in  particular  r.on- 
averaged  transects.  The  mean  vertical  temperature  profile 
is  characterized  by  a  mixed  layer  at  13°C  limited  by  a  sharp 
thermocline  (on  the  average  0.2°C/M).  On  the  shelf,  a  bottom 
mixed  layer  of  1 2® C  was  present. 

Taking  an  arbitrary  "lsveL  of  no  motion”  at  150M, 
dynamic  heights  were  computed  for  aach  XBT  soundings  and  a 
horizontal  smoothed  contouring  was  nade  (Fig.  4.14).  The 
surface  geostrophic  speeds  corresponding  to  this  analysis 
are  cyclonic  around  a  stretched  ellipse  with  a  major  axis 
coinciding  with  the  continental  slope.  Northwest  currents 
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maximum,  vs.  a  returning  southeast  flow  of  .2S/S  eve-  » 
rise . 

A  vertical  modal  analysis  (Cf.  7I-B)  ,  yields  ‘'are- 
clinic  infernal  Rossby  radii  cf  deformation  of  4.5K!1  or.  -he 
shelf  and  25KS  in  the  deep  ocean.  (Cf  Vl-B) 

2.  Analysis  of  the  Variability 

The  spatial  variab-lity  of  the  thermocline  depth  was 
computed  using  a  subset  of  the  total  X3T  data,  which  was 
selected  to  avoid  unrepresentative  (sparce  and  dispersed) 
values.  The  variance  of  the  depth  of  a  given  isotherm  was 
calculated  in  each  rectangular  bin  of  1 0 KM  width.  Standard 
deviations  of  the  1  3°C  isotherm  depth  from  the  above  means 
were  computed  (Pig.  4.9)  .  The  strong  increase  in  the  values 
at  the  position  of  the  idealized  shelf  break  could  either 
mean  a  significant  increase  in  the  time  variability 
(infernal  waves)  ,  or  express  the  fart  that  the  mixing  is 
faking  place  over  a  very  short  distance  compared  with  th« 
width  (10KH)  of  the  bins  used  for  averaging. 

3.  Select ed  Represents  tiya  Transects 

Seme  specific  transects  were  then  considered  to  give 
a  quasi- instants neous  spatial  image  of  this  mixing  process. 
One  interesting  subset  of  BT’s  was  acquired  by  a  transiting 
ship  sailing  towards  Ouessar.t  (Osnant  Island)  and  whose 
track  crossed  the  shelfbreak  at  47<>30S,  06°W  on  August  23.  A 
sampling  rate  of  one  per  hour  (i.  =  .  30KM  spacing)  enables 
drawing  a  quasi-synoptic  transect  parallel  to  the  SAS  swath 
but  60KM  west  of  it  (Fig.  4.10).  Other  transects  were  made 
by  selecting  3T*s  close  in  rime  and  oonvenianfly  spaced  in 
the  vicinity  of  a  perpendicular  to  the  isobaths.  A  se*  of 
three  transects  is  depicted  on  a  geographic  grid  (Fig.  4.11) 
with  other  non-synoptic  transects  reproduced  from  various 
sources  (Fic.  4.12).  The  scaling  adopred  is  identical 


for 


all  and  bas*d  on  the  "alon g-*rack"  distance  from  the  shelf- 
break.  Common  characteristics  appear  zo  be  a  narrow  pelagic 
front  anchored  on  the  sea-floor  at  the  upper  part  of  the 
shelfbreak  (ca .  1803  depth)  ,  with  one  or  two  upwellirg  zones 
adjacent  to  this  narrow  intermediate  mixed  layer.  The  ther- 
mocline  is  stronger  on  the  shelf  than  in  deep  waters  (Fig. 
4.13).  In  chapter  711 ,  an  idealized  3-layered  ocean  was 
extrapolated  from  these  selected  transects  and  the  corre¬ 
sponding  geostrophic  currents  ware  compared  to  the  observed 
circulation . 

2.  TEMPORAL  THERMAL  ANALYSIS 

Daily  averages  of  the  isothermal  depths  (Fig.  4.15)  show 
a  significant  heating  in  the  upper  layers  on  20  August  ar.d  a 
general  broadening  of  the  thermocline  between  17  and  24 
August.  Assuming  a  propagation  offshore  from  the  shelf  of 
tidally  generated  waves,  diagrams  have  been  established 
displaying  for  a  given  isotherm,  variations  of  its  depth  as 
a  function  of  time  and  space  (Fig.  4.  16)  .  These  variations 
in  depth  were  obtained  by  subtracting  the  mean  along  the 
x-directicn  to  eliminate  any  permanent  trend  and  enable  a 
contouring  of  relative  values.  The  spatial  axis  is  th° 
along-track  (/-coordinate)  distance  from  the  shelfbreak;  the 
time  axis  in  hours  has  been  "folded  ov»r"  two  semidiurnal 
(12. 4H)  periods.  Hence  ail  time  coordinates  are  taken 
within  24  hcurs  by  subtraction  (or  addition)  to  the  actual 
time  of  observation  cf  an  adequate  multiple  of  cha  semidi¬ 
urnal  period. 

Mo  contouring  algorithm  could  oe  applied  to  such  sparse 
and  greatly  contaminated  data  (by  various  phenomena;  vari¬ 
able  wind  stress,  radiation  heat  exchang3,  ever,  high 
frequency  nonlinear  internal  waves).  A  hand-contouring  was 
attempted  under  the  following  orincioles:  one  contour  line 


was  drawn  for  a  zerc  relative  depth,  another  line  followed 
the  maxima  and  the  last  one  followed  the  minima.  This  was 
done  for  three  temperatures  (12,14 , 15®C)  .  Each  of  *hosa 
diagrams  exhibits  a  constant  pattern  obliquely  oriented  on 
the  axis,  and  thus  corresponding  to  a  propagation  of  the 
variations.  The  slope  of  those  patterns  (disnar.ce  of  propa¬ 
gation  during  a  given  interval  of  time}  was  found  to  be 
approximately  1.1M/S  (equivalent  to  D.  9M/S  perpendicular  to 
the  shelfbreak)  .  A  modal  analysis  described  in  VI-B  showed 
that  high  frequency  internal  waves  of  period  90MIN  would 
have  this  same  wave  speed.  Furthermore  the  phase  of  this 
propaganion  implied  a  maximum  deepening  event  at  "Y  =  0M  at 
about  1900Z,  20  Aug.  Conversely,  at  the  time  of  the  SAR 

image  the  crests  of  these  propagating  oscillations  were 
found  to  be  at  Y*10  and  Y=60.  Maximum  values  of  relative 
depth  excursions  show  general  increasing  values  with  depth: 

a.  For  isotherm  12°C  mean  dep*h  of  35M  with  55m  excursions 

b.  For  isotherm  14®C  mean  depth  of  45M  with  15M  excursions 

c.  For  isotherm  16°C  mean  depth  of  30M  with  17m  excursions 

At  deeper  depths,  the  NEADS  mooring  temperature  data 

allows  estimation  of  the  vertical  isotherm  movement  by 
relating  semidiurnal  temperature  fluctuations  at  a  given 
depth  to  the  climatological  temperature  profiles.  Data  from 
the  current  meter  pressure  gauge  showed  that  vertical  oscil¬ 
lations  of  the  array  itself  were  negligible. 

At  1500M,  temperature  variations  nf  .3®C  in  a  mean  gradient 
of  8  10'3  °C/M  yield  vertical  amplitudes  of  about  20M. 

At  3000M,  temperature  variations  of  . 15®C  in  a  mean  gradient 
of  2  10"*  ®C/M  yield  vertical  amplitudes  of  more  than  uom. 

At  3500M,  temperature  variations  are  at  the  limit  of 
discrimination  of  the  thermal  gauge  but  reveal  vertical 
movements  of  at  least  20M. 


F.  SUHHABT 

SEASAT  flaw  over  fha  shelfbreak  of  the  Bay  of  Biscay  at 
2142Z  (after  sunset)  or.  20  August  1373.  A  high  pressure  cell 
hai  dominated  the  area  since  early  in  the  morning,  brir.aina 
a  southward  wind  of  6M/S.  The  cloud  cover  gradually 
decreased,  and  the  air  temperature  decreased  to  17oc. 

A  well  established  swell  was  entering  th*  Bay  from  the 
northwest,  nearly  tangential  *o  the  isobaths  in  the  vicinity 
of  the  shelfbreak  where  a  significant  increase  in  wave 
heights  was  observed  from  the  SEASAT  altimeter. 

The  tide  was  a  spring  tide  at  high  level  plus  5  hours  at 
the  shelfbreak.  The  ebbing  surface  current  on  the  shelf 
side  was  at  its  maximum  toward  240°T  at  0.85M/S. 

The  ocean  was  strongly  stratified  with  sea  surface  temp¬ 
eratures  of  18°C ,  a  thsrmocline  gradient  of  3.3°C/?i  at  an 
average  dep*h  of  30n,  and  a  bottom  tidally  mixed  layer  on 
th9  shelf  of  12°C.  The  first  mode  internal  Rossby  radii  of 
deformation  ware  of  4.5KM  on  the  shelf  and  of  25km  in  deep 
water. 

A  spatial  thermal  analysis  showed  a  marked  front  at  the 
shelf  break  with  intense  mixing  and  upwelling.  There  was  a 
westward  surface  geostrophic  jet  at  the  shelfbreak. 

A  temporal  thermal  analysis  gave  some  evidence  of  an 
internal  semi-diurnal  "disturbance"  propagating  offshore 
from  the  break  with  a  spaed  of  3.3M/S  and  with  the  last 
crest  in  the  vicinity  of  "Y=1 0"  at  the  time  of  the  SAR 
image.  This  propagation  sp*ed  being  significantly  smaller 
than  the  computed  phase  speed  of  the  internal  first  modes 
tides  (3.8  and  1.7M/S),  either  some  higher  modes  are  presen* 
or  high  frequency  waves  are  generated  by  the  tide  but  even¬ 
tually  decoupled  from  if  and  trapped  in  the  vicinity  of  the 
thermocline  (where  the  largest  Brunt- 7a isala  frequencies  can 
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7.  SAR  IH&GE  ABALTSIS 

i.  OHBIT  P1BAMBTEHS  AND  GBOGRAPHIC  R  EFEBEHCES 

The  auxiliary  grid  introduced  previously  (1,4)  will  be 
used  to  conveniently  locate  features  with  respect  to  the  S&R 
image  and  to  the  average  shelfbreak  contour.  The  geographic 
fit  of  the  SAR  image  is  estimated  to  be  accurate  to  a*  least 
♦/-  5 KM  in  the  vicinity  of  the  shelf-break. 

At  2142Z  Aug  20,  SEASAT  was  on  a  descending  node  over 
Prance  with  the  SAR  antenna  oriented  to  -.he  ncrrhwest. 
Consequently,  a  positive  (northward)  shear  of  eastward 
current  would  appear  on  the  image  as  artificially  compressed 
with  enhanced  Bragg-scatt er  (Cf  Fig.  2.3)  In  oth*r  words,  a 
current  directed  toward  the  northwest  would  appear  on  the 
SAR  image  with  a  bright  northeastern  boundary  and  a  dark 
southwestern  one. 
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B.  HEURISTIC  ASALISIS  OP  THE  SAH  ISAGE 

The  patterns  visible  in  the  optically  processed  SAR 
image  may  be  grouped  in  four  distinct  categories  according 
to  their  geographic  location  and  to  their  geometric  charac¬ 
teristics.  They  are  discussed  separately  under  the  titles: 
"inner  shelf",  "outer  shelf",  "generation  area",  and  "propa¬ 
gation  area".  The  image  discussed  being  a  positive  print, 
the  lowest  baokscatter  radar  response  will  be  referred  *0  as 
dark  areas  and,  conversely,  the  enhanced  signal  returns  will 
be  referred  to  as  bright  or  light  areas. 

A  correlation  between  image  characteristics  and  bathyme¬ 
tric  features  is  made  and  must  take  inro  account  two  inaccu¬ 
racies.  First  the  error  in  mapping  the  SAR  image  is 
estimated  as  5KM  at  most.  SecondLy  the  bathymetric  chart 
used  here  (Fig.  5.5)  was  constructed  by  interpolating  sound¬ 
ings  of  different  sources,  with  gaps  of  ca.  10KM  in  some 
places  of  the  continental  slope,  dence,  the  relative  posi¬ 
tions  and  shapes  of  the  canyons  drawn  are  subject  to 
caution  within  the  same  order  of  magnitude  as  the  SAR 
feature  locations.  Nevertheless,  in  the  case  of  several 
assumed  correlations  between  well-defined  imaoe  and  bathyme¬ 
tric  features,  the  confidence  may  be  increased  by  the 
perfect  relative  fit  of  each  of  the  poinrs  (Fig.  5.4). 

1 .  Inner  Shelf 

(Y>40)  -  In  this  area,  manifestations  of  internal 
activity  are  of  low  contrast  and  lack  any  obvious  organisa¬ 
tion.  Separated  curved  lines  can  nevertheless  be  identified 
with  an  average  radius  of  curvature  of  12  K H.  Patterns  of 
similar  shapes  can  be  associated  an  spacings  of  28KM  (Fig. 
5.2),  but  -.heir  relation  to  any  kind  of  periodic  phenomena 
is  not  proven.  They  could  be  due  to  mesoscale  atmospheric 
variability,  or  linear  slicks  associated  winh  a  localized 
convergence  (If  para  II-C-2)  . 
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The  lowest  limits  of  wind  speeds  sufficient  fir  a 
radar  return  are  attained  in  the  northern  par*  of  this  area, 
particularly  south  of  the  "Chaussee  of  Sein"  and  leeward  of 
the  Britanny  Zoast  where  mesoscale  land-sea  breeze  may  well 
have  been  predominant.  (A  ship  sailing  100°r  can  clearly  be 
seen  with  its  wake  virtually  displaced  according  to  the 
Doppler  effect  described  in  II-B-1)  . 

2.  Out er  Shelf 

(10<Y<40)-  Internal  wave  group  patterns  are  apparent 
in  this  region.  They  appear  as  arcs  of  circles,  concave  to 
the  west  with  a  variable  radius  of  curvature  of  minimum 
value  15KM.  They  are  limited  to  the  north  by  a  caustic  in 
which  they  merge;  this  line  is  brighter,  i.e.,  of  greater 
backscatter.  They  are  broken  tc  the  south  by  another  line 
of  brighter  characteristics  still.  These  patterns  seem  to 
exist  over  the  entire  width  of  tae  3AR  image,  but  they 
appear  with  more  contrast  on  the  western  part  of  it.  Som= 
periodicity,  of  about  28KM  lengta  on  the  x- axis,  can  be 
detected. 

The  wavelength  of  these  patterns  is  variable  in  two 
directions.  (1) toward  the  so-called  caustic  to  the  north, 
the  wavelengths  decrease  to  zero,  starting  at  a  maximum 
value  on  the  southern  boundary.  (2)  toward  the  southeast 
(increasing  x-direction)  ,  which  is  their  apparent  direction 
of  propagation,  the  average  wavelengths,  measured  at  mid- 
distance  between  the  two  boundaries,  is  of  the  order  of  1KM 
on  the  western  side,  increasing  to  2 XX  in  th®  middle  of  the 
swath. 

The  most  conspicuous  wavecrests  are  located  at 
x=  - 2 5  ,  where  a  high  level  contrast  particularly  enhances 
those  curved  lines.  Each  of  thsse  arcs  appears  as  two 

narrow  contiguous  stripes:  a  dark  (low  backscatter)  one  on 
*he  concavity  and  a  bright  (high  backscatter)  one  on  th® 


convexity.  Hence,  over  one  wavelength  on  a  direc~ion 
perpendicular  to  the  crests  and  following  the  assumed  direc¬ 
tion  of  propagation  (from  the  inferior  to  the  exterior  of 
the  arcs),  the  backscatter  varies  approximately  in  the 
following  way:  7D%  of  the  length  at  a  uniform  level,  20%  at 
low  backscatter  level  and  10%  at  a  nigh  level  on  the  leading 
edge.  According  to  the  hydr odynamical  theory  of  SAR  imaging, 
this  result  means  that  the  internal  waves  responsible  for 
these  patterns  are  long-crested  with  deep  narrow  troughs 
(Pig.  2.4). 

The  southern  boundary  is  a  general  line  brighter 
than  the  nortaern  caustic  particularly  on  the  western  part 
of  the  image,  where  it  is  conspicuous  both  by  its  brighter 
characteristics  and  by  the  obvious  discontinuity  it  creates 
in  the  general  patterns.  This  boundary  is  located  approxi¬ 
mately  (♦/-  5KM)  over  the  180M  isobath,  thus  corresponding 
to  the  upper  break  in  the  continental  slope. 

3.  Contin ental  Slope  2*  Generat ion  Area 

(0< Y<1 0)  -  This  area  covers  the  upper  part  of  the 
slope  between  the  150  and  1000M  isobaths.  Once  again  the 
mesoscale  features  are  concentrated  oi  the  western  side  of 
the  image.  The  average  backscatter  is  weaker  (darker 
shading)  than  on  both  adjacent  areas  (shelf  and  deep  ocean)  . 
The  northern  limit  is  the  band  of  lighter  characteristics 
described  above. 

The  limit  to  the  south  is  (conversely)  a  discontin¬ 
uous  band  of  very  low  backscatter  abouo  1500M  wide,  which  is 
clearly  visible  from  x= -6  0  to  x=-2D  and  from  x  =- 5  to  x= 1 5 . 
The  straightness  of  this  line  (vs.  the  other  oatterns)  and 
the  fact  that  if  is  not  bordered  by  any  alternating  clearer 
soripes,  seem  to  indicate  that  it  is  due,  not  to  a  propa¬ 
gating  phenomenon,  but  to  a  stationary  (if  nor  steady) 
internal  feature.  These  could  be  svmpooms  of  a  surface 


kins  statical  theory  of 


SAP 


divergence  according  to  the 
imaging. 

In  this  area,  well-organised  coher=ni  circular 
patterns  suggest  cylindrical  spreading  of  internal  waves 
front  several  distinct  point  sources.  Assuming  that  the 
direction  of  propagation  of  a  given  wavepacket  can  be  esti¬ 
mated  both  by  the  curvature  of  the  crests  and  by  the  rela¬ 
tive  wavelength  dispersion  sense,  the  following  main  point 
souces  can  be  identified  and  given  names  for  further  refer¬ 
ence  . 


PI  centered  at  X=-42 


P2  centered  at  X=-2  5 


?3  centered  at  X=3 


?4  centered  afc  X=  4 5 


=  4  7  0  0  ON 
=335°  30W 

=  46052N 
=335  0  1 8W 

=  46°45N 
=005°02  W 

=  46  0  25  N 
=004O42M 


Point  PI  is  at  the  head  of  Penmarc'h  Canyon  at  900M 

depth . 

Point  P2  is  at  the  head  of  a  narrow  and  oblique  but 
deep,  unnamed  canyon. 

Point  P3  is  at  the  head  of  Suilvinec  Canyon  at  800M 

depth. 

Point  P4  is  at  the  head  of  3lavst  Canyon  at  900N 

depth . 

It  is  observed  that,  at  the  Odet  canyon  there  is  no 
evident  corresponding  internal  wave  pattern;  a  plausible 
reason  for  this  could  be  that  this  canyon  though  steep,  does 
net  go  up  to  the  shelfbreak  itself  (Fig.  5.5)  . 


At  pi  the  leading  wavs  of  the  nearest  group  is  15KM 
from  the  center.  Wavelengths  from  the  leading  trough 
decrease  from  2KM  to  0.7KM  with  a  range  of  +/-  50".  Som® 
interference  can  be  seen  with  a  second  circular  pat-ern  of 
poorer  definition  centered  around  ?*1  { X=-  44 ,  Y=-  1 0)  . 

At  P2,  two  wave  groups  propagate  in  two  distinct 
directions.  The  first  group  propagates  toward  130°T  along 
the  X-axis  with  its  leading  wave  at  10KM  from  the  center. 
The  second  group  propagates  with  a  30°  inclination  offshore 
(160°T)  ,  its  leading  wave  being  at  13KK  from  the  center  and 
having  the  greater  wavelength  in  the  group  as  above. 

At  P3,  no  obvious  group  can  be  identified  bu4-  the 
pattern  appears  as  a  superposition  of  faint  ellipses  with 
their  major  axis  oriented  130°T  along  the  X-axis. 

In  the  vicinity  of  Ph,  patterns  are  still  more 
confused,  but  a  guasi-circular  group  extends  around  if  at  a 
distance  of  17KH. 

h.  Continental  Rise  or  Propagation  Area 

(-90<r<0)  The  surface  patterns  here  are  more  blurred 
than  in  the  three  areas  described  above.  Features  are  on  a 
larger  scale  and  do  not  present  the  same  east-west  vari¬ 
ability.  The  brightness  is  uniform  over  the  whole  area.  The 
general  geometric  organization  is  that  of  arcs  of  circles 
roughly  centered  on  the  point  sources  defined  above,  or  more 
generally,  on  th»  1000M  isobath.  Some  alternating  black  and 
white  bands  may  be  identified,  but  are  larger  and  do  not 
contrast  sharply  on  a  grey  backdrop  (in  contrast  to  the 
so-called  generation  area).  The  distribution  of  those 
circular  bands  is  not  always  clear,  but  some  particular 
patterns  with  similar  shapes  can  be  found  at  a  spacing  of 
ca.  23  KM. 
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C.  AH&LYSIS  3P  THE  SBELL 

The  remotely  generated  swell  is  long  and  well  organized, 
thus  giving  generally  well  identified  patterns  on  th®  SAP. 
image.  Some  apparent  geographic  variability  in  direction  and 
crest  length  can  be  observed  and  could  be  indicative  either 
cf  variations  in  the  nonlinear  imaging  mechanism,  or  of  real 
deformations  in  the  swell  field  by  interaction  with  surface 
currents  or  mesoscale  winds.  The  wavelength  of  the  swell  is 
237M  (at  U7on,  004° .  30*  W)  ♦/-  30M.  This  corresponds  1 0  a 
period  of  12.3  seconds,  a  phase  speed  of  19.2  M/S ,  and  a 
group  speed  of  9.5  M/S.  The  direction  of  propagation  is 
well  defined  and  easily  measurable  on  the  shelf  and  offshor® 
Spain  on  the  southern  part  of  the  SAP  image.  In  the  vicinity 
cf  the  Armorican  shelfbreak,  the  patterns  are  sore  confused 
and  mix  with  the  internal  wave  field.  On  the  shelf,  the 
general  orientation  is  115<>t  (at  47<>if,  004<>W)  .  in  the  open 
ocean  (46°N,006°W)  ,  it  is  1350T.  South  of  the  shelfbreak 
if  seems  to  drastically  bend  towards  the  south  as  if 
refracted  by  an  opposing  current.  This  particular  feature 
will  be  discussed  further  in  section  7T-8. 

D.  SUMMARY 

On  the  shelf,  and,  particularly  in  the  first  30  KM 
inshore  of  ■‘■he  shelfbreak,  the  general  orientation  of  the 
surface  patterns  suggest  an  internal  motion  propagating 
along  the  isobaths  and  foward  the  east.  A  strong  and  local¬ 
ized  divergence  follows  ‘he  shelfbreak  at  180M  of  deDth, 
clearly  limiting  to  the  north  a  generation  area  on  the  slope 
where  several  point-sourcec  can  be  identified  at  *h®  heads 
of  sharp  canyons.  Internal  waves  are  generated  in  packets 
which  are  clearly  visible  at  short  distances  from  *:h®ir 
generation  points.  Wavelengths  of  the  internal  waves,  as 
well  as  the  spacing  between  packets,  are  highly  variable. 
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▼I.  CORRELATIVE  INTER PBETATI3B  OP  THE  SAR  IMAGE 


The  entire  area  of  the  continental  slope  imacrad  by  SAR 
is  dominated  by  surface  manifestations  of  short  internal 
waves.  These  internal  waves  appear  in  groups  which  are  prob¬ 
ably  generated  at  the  shelfbreak  in  association  with  a 
semidiurnal  internal  tide,  but  eventually  decouple  from  if 
as  both  perturbations  propagate  offshore.  Two  different 
fields  of  infernal  waves  have  to  be  considered;  (lithe 
internal  tidal  wave  with  a  semidiurnal  period  of  12.4  hours, 
and  (2)  the  higher  frequency  waves  generated  in  packers  as 
forced  by  the  internal  tide  interacting  with  shelf  bathy¬ 
metry,  or  by  dissipating  energy  at  higher  wavenumber  by 
breaking  processes.  The  characteristics  of  these  internal 
waves  are  derived  from  their  surface  manifestations  (Cf 
chapter  V).  Zomparing  the  results  with  a  synoptic  conven¬ 
tional  analysis  (Cf  chapter  17)  of  thermal  structure  and 
Eulerian  current  measurements,  gives  a  notion  of  the 
validity  of  an  operational  approaci  where  a  SAR  image  could 
be  the  basis  for  hindcasting  the  actual  internal  motions. 

A.  TWO  LAYER  MODEL  ON  THE  SHELP 

The  sharpness  of  *he  thermocline  on  the  continental 
shelf  allows  the  assumption,  in  a  first  step,  of  a  two-layer 
ocean  and  use  of  the  corresponding  dispersion  relation  to 
compute  phase  speeds  for  given  observed  wavelengths.  The 
dispersion  equation  is; 

_9_  AP _ 
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Or.  the  outer  shelf,  for  a  giver,  internal  wave  crest,  the 
observed  increase  of  wavelength  from  the  "caustics"  (Cf 
V-B-2)  to  the  leading  middle  point  averages  25$.  Working  by 
logarithmic  differentiation  on  the  simplified  form  of  the 
dispersion  equation 


gives  an  approximate  relation  between  this  wavelength  rela¬ 
tive  variations  and  the  increase  in  mixed  layer  depth: 

2  -  2  ^  ^  ^  ^ 

c.  x  u  -  H  * 

and 

{ jk .  2. 1L  ,  . 

X  U  -  7.  K 

According  to  the  observed  variations  ar.d  to  the  mean 
thermocline  depth  derived  from  XBT  data,  the  thermccline 
depth  is  estimated  to  vary  or.  the  outer  shelf  from  25M  on 
the  boundaries  to  55M  at  mid-distance  from  each  of  th^rn. 

On  the  continental  slope,  a  linear  interfacial  wave  of 
wavelength  1  KM  yields  a  phase  speed  of  0.48«/s  whereas  for 
700s  tne  phase  speed  decreases  to  0.47S/S.  Applying  these 
results  to  the  distance  separating  the  wave  packers  gener¬ 
ated  on  the  shelf,  and  assuming  tie  generation  points  are 
those  bathymetric  features  described  above  (Cf.III-B),  the 
travel  time  would  be  9  hours  for  the  leading  waves,  thus 
implying  a  generation  at  low  tide  with  maximum  ebbing 
surface  currents  (Cf  Table  4.6). 

Conversely,  and  if  the  generation  time  was  known  with 
better  precision  relative  to  the  tide  phase,  observation  of 
the  position  of  similar  groups  together  with  density 
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gradient  data  could  give  an  estimate  of  the  *■  har.no  cl  in  a 
depth . 

This  two-layer  approximation  would  hold  if  the  thermo- 
cline  was  sharp  enough  or  equivalently  ifr  at  its  level,  *he 
Brunt-Vaisala  frequency  profile  exhibited  a  strong  peak, 
thus  trapping  the  energy  of  high  frequency  internal  waves  in 
an  horizontal  channel.  But  indeed,  the  actual  ocean  has  not 
such  a  simplified  structure  ,  particularly  in  the  "genera¬ 
tion  area",  where  dramatic  mixing  processes  are  talcing 
place.  Furthermore,  ^at  low  frequencies  such  as  semidiurnal, 
the  vertical  movements  extend  to  tie  whole  water  column  and 
cannct  be  described  without  a  vertical  modal  analysis. 

B.  HOBHAL  BODES  OF  THE  BAROCLINIC  INTERNAL  TIDE 

Internal  baroclinic  modes  of  the  ocean  were  computed 
from  the  density  profiles  extracted  from  observed  temoera- 
ture  soundings  and  climatological  salinity  data  (see  Table 
71 II)  .  The  Br unt-vaisala  frequency  is  the  free  oscillation 
frequency  for  a  given  density  gradient  : 

n2=  - _2_(iLL_i_)  . 

$  '  Az  cL  • 

Density  gradients  were  calculated  by  the  equation  of  state 
from  the  observed  temperature  profiles  sampled  by  XBT  data, 
and  assuming  a  constant  vertical  salinity  profile  of  35.6%o 
in  the  deeD  water  and  of  35.<*%c  on  the  continental  slope 
from  climatological  data.  Brunt-7aisala  frequencies  were 
then  computed  in  the  first  300B  and  eventually  blended  to  a 
deep  mean  climatological  profile  found  in  Gould  ’  1973], 

The  eigenfunctions  solving  the  Sturm  Liouville  problem 
of  vertical  modes  give  the  vertical  dependence  of  the 
velocity  fields  for  any  frequency  of  internal-inertial 
oscillation.  An  infinite  number  of  solutions  exist  yielding 
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vertical  distributions  characterized  by  one  or  more  phase 
reversals  in  horizontal  speeds  with  depths  (Fig  6.1  and 
6.2).  For  eacn  of  these  modes  internal  wavelengths  and  phase 
speeds  can  be  found  as  a  function  of  the  fcroing  frequency. 

Computations  were  accomplished  using  a  numerical  solu¬ 
tion  [Mooers  1971]  to  the  eigenvalue  problem.  Starting  with 
the  boundary  condition  a*  the  surface,  successive  correc¬ 
tions  of  *he  eigenvalue  are  made  at  each  iteration  to 
satisfy  the  bottom  boundary  condition. 

The  solutions  are  defined  within  a  multiplicative 
factor.  An  estimate  of  the  vertical  distribution  of  ampli¬ 
tudes  values  is  then  necessary  to  resolve  this  incertitude. 
A  simple  best  fit  of  the  amplitude  values  from  available 
data  was  made  to  the  first  baroclinic  mode.  Vertical  ampli¬ 
tudes  at  the  semidiurnal  frequency  were  extracted  from  the 
current  meter  data  (Cf  IV-C)  :  temperature  fluctuations  as 
related  to  the  mean  vertical  thermal  gradient  at  the  corre¬ 
sponding  depth  gave  an  estimate  of  the  vertical  excursions 
of  the  isotherms  (as  shown  in  IV-Ei  .  (The  variations  of  the 
current  meter  pressure  gauge  showed  that  vertical  oscilla¬ 
tions  of  thc  array  itself  were  negligible) 

The  second  mode  vertical  profile  is  coherent  with  those 
observations  and  is  consequently  assumed  to  be  the  mean 
component  of  infernal  energy  at  the  semi-diurnal  frequency 
(Fig.  6.3).  Maximum  vertical  amplitudes  are  inferred  to  be 
at  2000M  and  of  the  order  of  50M,  a  level  of  no  vertical 
movement  (medal  node)  exist ina  at  1333M. 

The  analysis  of  the  BT  soundings  exhibits  an  offshore 
spacing  of  perturbations  of  45KM  corresponding  fc  phase 
speeds  of  1M/3  for  a  semidiurnal  generation  and  wavelengths 
of  7.2KM  in  the  open  sea  (90MIN  period!  from  the  first  mode 
dispersion  analysis. 


at  least  in  its 


Ir.  conclusion,  the  internal  tile, 
mole,  does  not  appear  an  the  SAR  image,  neither 
spatial  analysis  of  the  thermostructure.  Nevertheless 


odic  disturbances  could  be  generated  at  a  semi-diurnal  rate. 


consisting  of  high  frequency  waves  which  are  discussed  now. 


TABLE  VIII 
INTERNAL  WAVELENGTHS 


WAVELENGTH  (KM) 


—  -1 

1 

-1 

PERIODS 
(Hours)  I 

M0  DEO  1 
tropic) 

M0 DEI  i 

73. 

\  MODE 2  ' 
Baroc] 

M0DE3  1 
.inic  me 

|  MODE  4 
5des 

SHALLOW  WATER 

1  50M 

12.4  | 

6.0 
1.2  1 
0.4  i 

2570.  j 
8  88  . 
166.  | 

27.3  3 
9.30 
1.73 
0.53 

10.00 
3.50  ! 
3.63  1 
0.19  ; 

5.60 

1  .90 
0. 35 
0.  12 

3.50 

1  .20 
0.22 
0.08 

MIXING  AREA 
220M 

12.4 
6.0 
1.2  1 
0.4 

-  | 

3110.  | 
10  76. 
201.  | 
67.  | 

L  j 
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C.  HIGH  PREQOBNCY  INTERNAL  SAVES 


The  wavelengths  in  the  internal  wave  packets  are  of  the 
oriar  of  3 KM,  yielding  frequencies  of  the  order  of  70 HI 5 
which  are  consistent  with  the  internal  wave  climatology  (Cf 
III-  F) .  The  spacing  of  the  high  frequency  wav®  packers  as 
seen  in  the  propagation  area  of  the  SAR  image  is  approxi¬ 
mately  23KM. 

On  the  continental  shelf,  *  he  periodic  spacing  in  the 
patterns  is  about  28KM,  which  is  olose  to  the  semidiurnal 
internal  wavelength  computed  for  shallow  water. 

1 .  Estimation  of  the  Mixed  Layer  Oeoth 

It  is  assumed  that  the  climatological  density  struc¬ 
ture  of  the  Bay  of  Biscay  is  known,  and  that  frequent  XBT 
surveys  as  well  as  remote  sensing  of  the  sea  surface  temper¬ 
ature,  allow  estimation  at  any  time  of  the  mean  density  jump 
a*  the  thermocline.  Starting  vita  only  this  hypothetical 
element  and  making  abstraction  of  the  actual  in-situ  data 
available,  it  is  attempted  to  compute  the  thermocline  depth 
as  a  function  of  the  geomaf  t  characteristics  of  internal 
waves  as  derived  from  the  SAR  image. 

Assuming  a  climatological  density  jump  of  1.5x10  ’ 
at  the  thermocline,  and  observing  separation  between  wave 
packets  of  23KM  in  the  open  sea  area  of  the  SAR  image  yields 
an  estimate  of  the  mixed  layer  depth  of  20M.  A  variation  of 
10%  in  the  assumed  density  jump  would  result  in  10%  differ¬ 
ences  in  this  estimate,  whereas  changes  of  10%  in  the  inter¬ 
packet  distances  would  result  in  23%  changes  in  this  derived 
mixed  layer  depth.  Consequently,  if  the  absolute  estimate 
of  the  thermocline  d9pth  is  uncertain,  comparison  between 
interpacket  distances  would  be  a  very  sensitive  way  of  esti¬ 
mating  time  variations. 


va  ves 


2.  Estimation  of  the  I  nternal  »a  ve_s  Am  pi  it  ude 


Amplitudes  of  the  nonlinear  short  internal  waves 
wars  computed  by  a  method  described  and  applied  on  a  similar 
SAR  image  by  Apel  and  Gonzalez  [  1983  ].  Tha  first  assumption 
made  was  that  in  the  "generation  area"  (continental  slope), 
tha  internal  waves  visible  on  the  SAR  image  could  reasonably 
be  considered  as  solitons.  By  definition,  a  soliton  is  a 
perturbation  oonsarving  its  shape  and  general  characteris¬ 
tics  (excep"  phase) ,  even  after  collision  wi^h  another  one. 
This  characteristic  is  apparent  in  the  vicinity  of  point  Pi 
for  example  (Pig.  5.3)  ,  where  several  wave  packers  are 
superimposed  without  any  notable  mutual  interference. 

This  conservative  nature  of  solitons  is  linked  to 
tha  fact  that  the  linear  phase  speed  dispersion  is  counter¬ 
acted  by  a  nonlinear  acceleration  dependent  on  the  amplitude 
itself.  This  balance  is  expressed  in  the  last  two  terms  of 
the  Kortweg  -  de  Tries  equation: 


"4  +  , 

where  is  the  amplitude  of  the  solitary  wave  propagating 
with  time  "t"  in  the  "x"  direction,  "Y"  is  the  coefficient 
representing  the  linear  dispersion  and  "«C"  is  the  nonlinear 
coefficient  and  is  given  by 


where  " W "  is  the  vertical  eigenfunction  or  structure  func¬ 
tion,  "h  ,  "  and  "h^"  ice  the  depths  of  the  layers  in  a 
simplified  discrete  model. 

From  the  mean  density  stratification  observed  at  the 
time  of  the  SAR  image,  and  using  the  vertical  modal  analysis 
described  above,  linear  phase  speads  "C  "  for  each  observed 
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wavelength,  and  a  phase  speed  limit  "C  ”  for  infinitely  leng 
internal  waves  in  a  non  rotating  ocean  (steady- state  solu¬ 
tion  to  the  Sturm-Liouville  problem  were  compared. 

In  one  given  packet  ,  the  wavelength  stretch,  ("(")# 
was  defined  as  the  ratio  of  the  observed  wavelength  no  the 
linear  theoretical  one.  From  the  theory  [Apel  and  Gonzales 
1993],  this  ratio  is  equal  tc  twice  the  complete 

elliptic  integral  of  the  first  kind  with  argument  m2;  the 
argument  m2  bainc  itself  a  function  of  the  solitcn's  ampli¬ 
tude  A: 

(w1 --£*=£- A. 

C6-CP  3 

For  any  wavelength  it  is  then  possible  to  extract 
the  nonlinearity  due  to  a  given  amplitude.  Dr,  conversely, 
observing  a  nonlinear  dispersion  in  the  patterns  of  the 
internal  waves  (by  measuring  the  variation  of  wavelength  in 
a  given  packet) ,  it  is  possible  to  derive  an  ®srimated 
amplitude  of  the  waves. 

was  taken  to  be  of  tha  order  of  5x10*2  m  *,  ("<*.”  "sIOM)  . 
" C0"  was  found  to  be  0.55M/S. 

For  each  observed  wavelength  in  a  given  packet,  the  corre¬ 
sponding  phase  speed  can  be  computed,  and  the  wavelength 
stretch  can  be  measured.  A  value  for  "m2"  is  derived  from 
which  wave  amplitude  is  calculated  as 

A  -  c<>- g  P  3>  m1 

For  an  average  stretching,  taken  from  the  wavegroups 
visible  in  rhe  generation  ar®a,  amplitudes  of  15M  for  th® 
leading  wav«s  were  found. 
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D.  INTERNAL  IA7E  INTERACTION  WITH  5BELFBREAK  STOCTORB 

Two  different  kinds  of  geographic  changes  in  the 
internal  wave  group  characteristics  are  apparent  (Cf  7-3-3). 
The  general  shortening  in  wavelength  will  be  assumed  here  to 
be  either  a  Doppler  shift  due  to  a  ourrer.t  shear  or  a  phase 
speed  compression  due  to  a  horizontal  density  gradient  (e.g. 
as  in  a  front  or  a  mixing  area) .  Refraction  of  the  crests 
resulting  in  a  deformation  from  their  original  cylindrical 
spreading  as  in  groups  northeast  of  Pi  ar.I  ?2,  can,  on  the 
ether  hand,  be  more  surely  associated  with  a  strong  near 
surface  current.  Estimates  of  the  direction  of  a  shear 
surface  current  as  given  by  characterise!:;  refraction 
patterns  are  indicated  on  figure  5.4  They  give  some  indica¬ 
tion  of  a  general  southwest  ward  alongshore  flow  to  the  south 
of  the  generation  area  (on  the  continental  slooe)  ,  opposite 
to  the  more  localized  currents  at  the  shelf  break,  and  to 
those  south  of  the  slops  (Pig.  6.4*. 

E.  SWELL  REFRACTION  AND  SBELFBREAK  CURRENT 

A  subjective  analysis  of  the  swell  refraction  was 
attempted.  No  Fourier  Transforms  were  made  of  the  SAR  image, 
but  in  most  parts  of  if,  the  clarity  of  the  swell  field  is 
such  that  visual  estimation  of  its  direction  and  (with  less 
precision)  cf  its  length  can  be  made.  It  was  thought  inter¬ 
esting  to  try  to  apply  recent  theories  on  wave/current 
interaction  to  this  particular  case,  where  some  severe 
refraction  is  obvious  and  the  presence  of  an  interacting  jet 
probable.  The  theory  initially  developed  by  Phillips 
[1981],  has  been  applied  to  a  SAR  image  of  the  Gulf  Stream 
by  Seadows  [  1983  ];  it  was  transposed  with  drastic  simplifi¬ 
cations  to  this  case. 


In  the  sane  way  "■-hat  the  Snell’s  law  (lot  de  Descartes), 
describes  the  refraction  of  a  ray  encountering  a  gradient  of 
propagation  in  a  medium,  the  following  equation  describe? 
the  refraction  of  swell  of  a  given  period,  propagating  with 
a  phase  velocity  "c"  through  an  interface  between  two  sea 
surface  areas  with  surface  currents  of  "uQ"  and  "U,"; 
and  "K\"  are  the  angle  of  incidence  and  of  refraction  of  a 
swell  "rayM  (Pig.  6.5). 
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Assuming  that  the  cause  of  the  observed  refraction  is  a 
shear  current  directed  along  the  shalfbreak  (or  at  least 
parallel  to  the  isobaths)  ,  the  angle  of  incidence  and  the 
angle  of  refraction  (just  after  the  first  visible  break  in 
the  wave  crests)  were  measured.  Tie  current  speed  (equiva¬ 
lently  the  jet  velocity)  was  estimated  to  be  0.6M/S.  The 
accuracy  obtained  cannct  be  more  than  50%,  but  the  basic 
observation  sesms  reliable  enough  ,  at  least,  to  suggest  the 
presence  of  a  jet  and  to  localize  the  position  of  the 
maximum  speed  at  the  continental  break. 


fix.  TBNT  ATIVS  EXTRAPOLATION  JO  3EHBRAL  OCEANIC  FEATURES 


A.  INTERNAL  ff  ATE  GENERATION 

"In  the  coastal .  region  internal  tides  generated  on  a 
continental  slope  are  standing  waves"  [Rattray  1959  ].  The 
modal  analysis  relative  to  the  density  structure  observed  at 
the  time  of  the  SAR  image  shows  wavsleng-hs  of  the  order  of 
27km  for  the  semidiurnal  first  mode  internal  tide  on  the 
shelf.  This  wavelength  is  comparable  with  the  spacino  of 
patterns  visible  on  the  SAR  image  north  of  fche  shelfbrealc. 
It  can  therefore  be  assumed  that  these  patterns  correspond 
to  the  interaction  of  another  internal  perturbation  periodi¬ 
cally  modulated  in  phase  with  ridges  and  troughs  in  an 
internal  tide.  Obviously,  the  assumption  of  a  perfect 
standing  internal  tide  would  only  be  true  if  the  boundaries 
of  the  shelf  3ea  acted  as  perfect  reflectors. 

Wavelengths  of  the  first  internal  baroolinic  modes  in 
deep  wat«r  ace  170KM  for  mode  one,  76KM  for  mode  two  and 
44 KM  for  mode  three.  No  direct  manifestation  at  these 
length  scales  is  visible  on  the  SAR  image.  Nevertheless,  the 
high  frequency  internal  wave  packets  give  an  indirect  demon¬ 
stration  of  the  high  energy  loss  encountered  by  the  internal 
tide  at  the  irregularities  in  the  shelfbrsak.  Located  at  the 
heads  of  the  narrow  canyons,  oscillations  in  the  water 
column  are  estimated  (for  this  specific  case  of  the  spring 
tide)  to  be  of  the  order  of  15H  which  are  assumed  *o  be  lee 
waves  forced  by  the  ebbing  tidal  current.  The  internal  wave 
packets  generated  by  this  mechanism  are  localized  at  the 
canyons  extending  toward  the  top  of  slope  tc  a*  least  the 
1090M  isobath  and  occur  in  bursts.  They  eventually  develop 
arisotropically  due  to  the  presence  of  a  shelfbreak  iron4-  to 
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the  north.  The  general  direction  of  propagation 
the  south-west  perpendicular  to  the  continental  slope.  A 
thorough  discussion  of  ths  possible  generation  mechanisms  is 
given  by  Pingree  ‘1981].  The  specific  influence  of  narrow 
canyons  has  been  recently  enlightened  by  a  laboratory  exDer- 
iment  [Baines  1983]. 

B.  TIDAL  Hill NG  OH  THE  COHTIHEMTAL  SLOPE 

Tidal  mixing  on  the  continental  slope  is  demonstrated  by 
the  divergence  in  isotherms  in  the  quasi-synoptic  XBT  anal¬ 
ysis,  (Fig.  ^ .  1 0  and  tt.13)  .  and  is  in  good  agreement  with 
other  measurements  made  in  summer  and  fall  at  the  same  loca¬ 
tion.  In  Pingree  (1981),  the  various  hypothesis  related  to 
this  mixing  are  developed. 


C.  SHELFBREAK  PROHTS  AND  JETS 

A  new  insight  on  the  shelfbrealc  circulation  is  given  by 
the  SAR  image.  Parallel  bands  of  divergence-convergence 
which  compress  internal  waves,  can  be  associated  with  the 
thermal  structure  inferred  from  BT  data  and  illustrated  by 
the  following  simplified  scheme  of  circulation  (Fig.  7.1). 
Over  the  shelfbrealc,  the  dynamic  heights  fend  to  be  lower 
than  on  either  side  as  a  consequence  of  the  tidal  mixing. 
Therefore,  a  strong  geostrophic  flow  extends  to  the  thermo- 
cline,  follows  the  isobaths,  and  is  confined  to  the  shallow 
part  of  the  continental  slope.  las  flow  is  consequently 
toward  the  northwest  on  the  upper  part  of  the  slope  and 
toward  the  southeast  in  the  offshore  half.  The  idealized 
transect  (Fig  7.1),  accounts  for  the  localized  and  strong 
northwest  flow  and  explains,  by  stretching  of  the  thcrmo- 
cline  offshore,  the  lower  values  estimated  for  the  counter¬ 
current.  It  is  nevertheless  obvious  that  any  kind  of 
analytical  model  based  on  geostrophy  (steady-state)  and 
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alongshore  homogeneity  rannct  be  fully  appropriate  whet 
applied  to  *he  Celtic  Sea  shelfbreaic.  It  this  area  current® 
are  dominant  and  limited  geographically  by  the  coast  of 
Britanr.y,  bathymetric  irregularities  are  essential  in  the 
internal  generation  processes,  and  nonlinear  effects  between 
the  steady-state  circulation  and  the  tidal  movements  may  be 
expected. 

The  most  informative  indication  of  the  presence  of  a 
strong  alongshore  jet  are  the  swell  refraction  and  the 
sudden  increase  in  wave  height  (IV-B,  7-C,  71-5)  .  The  coin¬ 
cidence  of  these  phenomena  give  reliable  evidence  of  a  jet 
which  may  be  related  to  the  poleward  eastern  boundary 
current  discussed  by  Swallow  [  1977]. 


106 


VIII.  APPLICATION  TO  OMDERJATBR  ACOgSTICS 


Of  the  host  of  variables  affecting  sound  propagation  in 
the  ocean,  the  temperature  profile  has  the  strongest  influ¬ 
ence.  Not  only  the  vertical  gradient  profiles,  but  also  the 
horizontal  gradient  fields,  determine  the  widely  varyino 
shapes  of  •‘the  insoni  fication  volumes. 

Hence  temporal  and  spatial  temperature  variability 
should  be  taken  into  con sideratinn  when  forecasting  any 
acoustic  quantity,  from  one  single  climatic  or  observed 
‘■emperature  profile.  Two  approaones  are  possible  in  “his 
matter:  (1)a  deterministic  approach  would  be  possible  only 

in  the  cases  where  the  water-mass's  vertical  movements  are 
known  as  functions  of  time  and  geographical  locations,  well 
enough  to  allow  the  forecast  of  effective  variations  of  the 
acoustic  transmission  loss  with  time,  location  and  orienta¬ 
tion  of  the  sound  beam,  ar.d  (2)  a  probabilistic  approach 
would  attempt  to  model  the  oscillations  without  any  phase- 
lock  and  would  only  yield  estimates  of  the  acoustic  propaga¬ 
tion  variability. 

Obviously,  in  the  common  case  where  the  geophysical 
variables  vary  on  a  short  time  scale  and  without  any  period¬ 
icity,  the  probabilistic  approach  is  the  only  one  possible. 
Or.  the  other  hand,  with  larger  time  and  space  scales,  an 
effort  in  sampling  the  actual  geophysical  environment  is 
eminently  beneficial.  In  the  SuLf  Stream  region,  for 
example,  knowledge  of  the  positions  of  cold  and  warm  eddies 
has  proven  instrumental  for  acoustic  forecasts.  In  the  case 
of  shorter  scales,  but  where  periodic  motions  are  related  to 
a  known  generation  mechanism  (e.g.,  tides),  the  same  kind  of 
hindcasting  may  be  possible. 
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It  is  a~t=mptad  here  to  show  how  the  ocean  variability 
observed  ir.  the  Bay  of  Biscay  ma>  affect  sonar  detect ior. 
ranges,  and  to  assess  a  possibility  of  analysing  the  depen¬ 
dence  of  these  ranges  on  ralativa  positions  of  source  and 
receiver  vs.  the  internal  propagating  wave  field,  as  well  as 
on  the  direction  of  propagation  relative  to  the  internal 
wave  orientation.  The  internal  variability  taken  as  a  basic 
assumption  is  an  idealized  extrapolation  of  the  results 
gathered  from  the  study  of  the  SAX  image  on  August  20  1978. 
Therefore,  if  applies  only  to  this  giver,  particular  situ¬ 
ation,  and  with  a  wide  range  of  uncertainty. 

This  approach  may  prove  feasible  in  similar  cases  where 
internal  waves  can  be  occasionally  sampled,  either  from 
space  using  a  SAH,  or  by  adequate  in  situ  sampling  techni¬ 
ques  based  on  a  thorough  knowledge  of  fhs  internal  wave 
climatology. 

The  influence  of  an  internal  wave  field  on  acoustic 
propagation  is  especially  critical  at  high  acoustic  frequen¬ 
cies,  where  diffraction  and  leakage  are  minimal  and  where 
pure  geometric  ray  tracing  can  be  applied.  The  ir.sonified 
volume  will  then  drastically  vary  in  the  three  dimensions. 
Long  infernal  waves  will  act  as  long  scale  time  variations 
ir.  an  observed  or  predicted  SSP  (sound  speed  profile)  .  The 
effect  cf  short  waves  can  be  approached  stochastically  and 
lead  to  establishing  statistical  conclusions  on  the  range 
variability.  The  variation  in  propagation  will  be  dramati¬ 
cally  dependent  on  (1)fhe  relative  position  of  source  (or 
receiver)  with  respect  to  the  wave  field,  (2)  the  relative 
direction  of  sound  propagation  and  wave  field  orientation, 
ar.d  (3)  the  mean  SSP  and  its  variability  in  space. 

Two  cases  are  to  be  considered  .  ( 1)  A  long  (e.g.  tidal) 

internal  wave  field  for  low  frequencies  (long  scales)  to 
evaluate  variation  on  long  ranges  (of  ordsr  20 K )  cf  sound 
oro pagan icr. :  in  this  case  the  Paraoolic  Zquaticn  (PS)  lo  del 


has  been  used.  (2)  A  short  (as  seen  on  the  5  AR  image  or  the 
continental  si  ope)  free  internal  wave  field  at  high  frequen¬ 
cies  (short  scales)  to  evaluate  interaction  with  a  high 
frequency  acoustic  wave.  In  this  case  the  MTDJSA  model  was 


A.  IDEALIZED  WAVE  FIELDS 

The  acoustic  model  was  evaluated  assuming  that  in  a 
vertical  plane  the  sound  speed  profile  varied  from  a  mean 
profile  in  a  periodic  way  (Fig.  8.1).  This  assumed  profile 
was  drawn  from  ^he  XBT  data,  and  was  then  idealized  and 
smoothed  to  present  a  perfect  surface  mixed  layer,  a  sharp 
thermocline  and  a  climatological  deep  profile.  For  a 
constant  SSP  two  convergence  zones  are  apparent:  one  du*  to 
the  Mediterranean  outflow  at  1 100M,  and  another  due  *o  the 
pressure  increase  with  depth.  The  sharp  thermocline  was 
made  to  vary  in  depth  without  any  deformation  in  its  thick¬ 
ness  or  its  slope  (Fig.  8.1);  the  siape  of  this  dep-h  varia¬ 
tion  was  mcdeLed  as  a  periodic  curve  with  wavelengths  of  10 
and  2 OKM . 

Increasing  this  internal  wavelength  by  a  factor  "  k"  is 
equivalent  to  modeling  the  ray  tracing  in  a  plane  oriented 
with  respect  to  the  wavecrests,  at  an  angle  "A"  such  that 

sin  (A)  =1 /k  . 

It  was  then  possible  to  calculate  transects  of  the  three- 
dimensional  iasonif ication  (or  transmission  loss)  in  various 
relative  directions  from  the  source  and  draw  horizontal 
diagrams  a*  one  given  depth  (Fig.  3.3  and  8.4). 


B.  TBiSSHISSXOH  LOSS  DERIVED  PROM  THE  PE  MODEL 


The  parabolic  eguation  (PE)  nodal  developed  by  NO RD A 
[Brock  1978],  modified  by  NUSC  [Lae  and  Botseas  1982],  ar.d 
adapted  (with  acre  modif icationsi  for  use  at  the  Naval 
Postgraduate  School,  gives  a  normal  mode  solution  to  the 
wave  problem  in  the  form  of  ar.  elliptic  equation.  The  calcu¬ 
lations  made  with  this  model  must  be  considered  with  some 
caution  as  the  initialization  step  will  affect  the  near- 
fiald  results.  Principal  significance  should,  therefore,  be 
attached  to  PE  model  calculations  at  the  convarger.ee  zones 
(CZ)  ranges.  In  the  case  where  the  source  was  placed  over 
an  internal  wave  crest,  no  significant  directional  depen¬ 
dence  was  observed  and  the  transmission  loss  was  equivalent 
to  that  found  in  a  constant  SSP  environment.  Conversely,  in 
the  case  whera  tha  source  was  placed  above  an  internal  wave 
trough,  the  convergence  zones  were  found  to  be  enhanced  by  a 
varying  SSP.  Horizontal  diagrams  of  transmission  loss  were 
drawn  for  two  levels  for  the  receiver  (15  and  45M)  by  inter¬ 
polation  between  five  vertical  transmission  loss  curves  at 
various  directions  relative  to  the  internal  wave  field  (Fig. 
8.4)  . 

C.  RAT  TRACIS  6  WITH  MED  OS  A 

The  MEDUSA  Range  Dependent  Ray  Model  was  developed  by 
tha  Applied  Research  Laboratories  of  the  University  of  Texas 
at  Austin  under  NURDA  sponsorship  during  FY80.  MEDUSA  uses 
smooth  functions  (splines)  to  interpolate  between  the 
discrete  SSP's  (given  as  input)  and  solves  the  various  ray 
equations  numerically  for  a  limited  number  of  characteristic 
rays  (eigenrays)  .  Eigenrays  are  located  by  interpolation  in 
range  (not  in  depth  as  in  other  ray  tracing  models)  between 
rays  bracketing  the  receiver  in  range.  The  source  is 
assumed  to  lie  on  the  vertical  axis  of  a  cylindrical  cocrdi- 
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rats  system,  around  which  the  geophysical  variables 
symmetrical.  Sigenrays  are  located  by  int erpolatin g  be-wee.n 
rays  bracketing  the  receiver  in  range  (not  in  depth)  . 
(Hence  any  direct  modeling  of  horizontal  refraction  is 
excluded).  The  PROFTHP  program  permits  diagnosis  cf  profile 
artifacts  produced  in  the  geophysical  interpolation  routine. 
RATPLT  produces  ray  path  plots. 

This  ray  tracing  program  was  ased  in  the  same  way  as 
the  PE  model  to  compute  acoustic  transmission  in  vertical 
planes  at  various  angles  with  an  internal  wave  field  orien¬ 
tation.  But  here,  the  features  which  were  investigated  were 
at  shorter  ranges;  the  direction  dependence  exhibited  in  the 
horizontal  diagrams  (Fig.  8.3)  relates  to  an  eventual 
refraction  under  the  mixed  layer,  thus  cancelling  (for  rays 
travelling  at  an  angle  with  the  internal  wave  crests)  the 
"shadow  zone".  Results  are  only  purely  qualitative  and  would 
need  to  be  investigated  from  a  transmission  loss  point  of 
view. 

Medusa  was  used  also  to  input  the  horizontal  sound  speed 
profile  as  observed  in  the  same  idealized  internal  wave 
field  along  an  horizontal  in  the  direction  of  wave  propaga¬ 
tion  and  at  the  mean  depth  cf  the  thermocline  (Fig.  8.2). 
The  result  of  ray  tracing  in  the  norizontal  plane  will  be 
valid  only  given  the  assumption  that  che  rays  do  not  make 
any  depth  excursion  from  the  chooser  depth  of  the  horizontal 
SSP.  This  assumption  can  be  made  only  because  the  vertical 
ray  tracing  shows  that  there  is  a  surface  sound  speed 
channel,  sc  that  a  significant  part  of  the  acoustic  energy 
will  propagate  within  the  range  of  variation  of  the  mixed 
layer.  Above  the  the  thermocline  the  sound  speed  will  not 
exhibit  any  horizontal  gradient  and  no  refraction  will 
occur . 
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D.  GsHBBAL  RESULTS 

Transmission  losses  (TL)  in  an  idealized  internal  wave 
field  were  found  to  be  equivalent  to  the  TL  relative  to  a 
constant  SSP  environment  if  the  source  was  placed  over  a 
crest.  In  the  case  where  the  source  is  above  a  trough,  and 
for  internal  wavelengths  of  10 KM,  the  transmission  loss 
decreases  at  the  convergence  zone  by  10dB  for  an  insonifica- 
tion  perpendicular  to  the  crests  (Fig.  8.4)  This  anisotropy 
at  long  ranges,  even'  if  not  representative  of  any  real  ocean 
environment,  is  at  least  an  illustration  that  passive  detec¬ 
tion  in  a  long  internal  wave  field  will  be  dramatically 
direction-dependent  and  suggest  the  importance  of  taking 
this  factor  into  account  in  any  passive  detection  scheme. 

Ray  tracing  shows  that  there  is  a  possibility  of  insoni- 
ficatior.  of  the  "shadow  zone"  beneath  the  thermocline  with  a 
variable  SSP  structure.  This  would  affect  both  passive  and 
active  search  (or  escape)  tactical  patterns  afc  low  ranges 
(of  the  order  of  4000M) .  Ray  tracing  in  the  horizontal 
suggests  the  possibility  of  an  horizontal  duct  (particularly 
in  deep  and  narrow  internal  troughs*  ,  but  the  assumptions 
made  here  are  too  drastic  to  allow  any  conclusion  from  the 
elementary  run  made  with  Medusa.  Further  work  in  this  line 
should,  nonetheless,  help  to  tailor  use  of  any  elementary 
(single  SSP)  ray  trace  model  to  interpret  some  consequences 
of  a  variable  SSP  on  transmission  anisotropy. 
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IX.  §i3  is  IN  INTER  NAL  OCEAI  SAMPLING  SYSTEM 

The  SEASiT  experiment  was  the  first  specific  oceano¬ 
graphic  use  of  a  SARfrom  space.  Assessments  made  a  priori 
on  its  monitoring  capabilities  have  been  justified  by  the 
large  number  of  oceanic  internal  processes  it  observed.  In 
less  than  42  hours,  it  covered  an  area  of  103  million  KM2  of 
the  Earth’s  surface  (35  of  them  over  the  oceans),  [Pravdo 
1983 ]. 

The  main  specific  advantages  of  a  5AR  as  an  ocean  remote 
sensor  system  are  its  all-weather  capability  (in  particular 
with  respect  to  cloud  cover,  night  conditions  and  sea- 
state),  its  high  resolution  and  its  sensitivity  to  the 
hydrodynamical  surface  manifestations  of  internal  activity. 

The  obvious  difficulty  encountered  with  an  operational 
use  of  such  a  sampling  system  is  the  computation  time 
required  to  process  the  enormous  amount  of  data  (some  10* 
pixels  per  second)  .  In  the  case  of  an  orbiting  station  (as 
proven  by  the  5IR-A  experiment  on  the  1982  Shuttle 
missions)  ,  the  data  may  be  acquired  and  exploited  on  board 
and/or  kept  in  memory  for  a  further  ground-processing.  But, 
for  an  unmanned  satellite,  data  recording  and  exploitation 
require  that  an  earth  station  be  at  radio  receiving  distance 
from  it. 

A.  GB1EBAL  REQUIREMENTS  FOR  A  SAR  3CEAN  MONITORING  SYSTEM 

To  possess  this  particular  capability  to  sample  the 
internal  ocean,  a  5AR  system  must  have  some  specific 
characteristics. 


B.  II  SITTJ  CDHPLBHENTAR  Y  S INULTAHEOOS  PROBING 

To  analyse  characteristics  of  the  internal  features 
displayed  by  a  3AR,  the  mean  density  structure  bus*  be 
known.  It  has  been  shown  (Cf  chapter  VI)  that  the  mixed 
layer  depth,  the  amplitudes  of  high  frequency  internal 
waves,  and  the  location  of  fronts  can  be  inferred  from 
consequent  surface  patterns.  If  would  therefore  be  possible 
to  conduct  such  an  analysis  in  real  tine  if  some  X3T 
sampling  was  made  nearly  synoptically  and  could  be  blended 
with  climatological  data. 

Obviously  time  series  of  temperature  structure  would  be 
more  efficient  still  by  giving  synoptic  information  on  the 
frequency  spectrum  and  vertical  modes  of  the  internal 
motions.  A  reiote,  drifting  thermistor  array  as  described  in 
[N3SC  1980],  if  launched  in  an  appropriate  position  with 
respect  tc  the  a  priori  internal  wave  field,  could  give  both 
ar.  absolute  reference  to  radar  remote  sensing  results  and 
fill  in  the  gaps  between  two  consecutive  passages  of  the 
satellite.  (A  radar  transpondor  could  add,  at  low  costs,  a 
geographic  reference.) 

C.  OPERATIONAL  PROCESSING  AND  EXPLOITATION 

The  analysis  of  a  SAR  image  may  be  made  or.  two  different 
levels  corresponding  to  different  goals,  time  constraints 
and  processing  facilities  (Fia.  9. It. 

A  heuristic,  manual  analysis  , similar  to  the  interpreta¬ 
tion  by  meteorologists  of  a  visible  or  infra-red  satellite 
image  of  the  atmosphere,  would  probably  be  effective  if 
other  archived  images  were  available.  A  simple  comparison  of 
similar  features  is  capaale  of  giving  indications  of  the 
relative  properties  of  a  given  internal  wave  field,  and  of 
the  geographic  displacement  of  fronts.  Rough  quantitative 
measurements  of  interpacket  separation,  wavelengths,  and 
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variation  of  wavelength,  could  b®  used  in  simple  ar.aly-icai 
models  to  compute  features  of  operational  interest  such  as 
internal  wave  amplitudes  and  directions,  and  location  of 
current  jets. 

On  the  other  hand,  for  routine  exploitation  of  5AR 
images  of  a  given  area,  the  knowledge  of  climatology  could 
probably  help  to  tailor  processing  algorithms,  eventually 
using  the  raw  data  directly.  Once  the  features  of  interest 
are  defined  by  their  approximate  locations,  dimensions  and 
average  characteristics,  the  reguired  data  could  be 
extracted  as  a  subset  and  processed  according  to  specific 
methods.  For  example,  if  a  front  is  known  no  move  parallel 
to  itself  and  to  be  characterized  on  a  SAR  image  by  a  band 
of  high  radar  intensity  due  to  convergent  enhancment  of  the 
surface  waves,  the  mean  radar  intensity  averaged  parallel  to 
the  front  would  be  enough  to  locate  its  position  (and  inten¬ 
sity)  .  If  the  presence  of  high  fraguency  internal  waves  is 
assumed,  a  procedure  to  extract  their  mean  direction  and  the 
distribution  of  their  wavelengths  would  be  to  analyse  a  set 
of  pixels  over  a  limited  area  by  a  two-dimensional  Fast 
Fourier  Transform  technique. 
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X.  SUMMARY  AND  RBCOagENDATIONS 


The  SAR  linage  gathered  by  SEASAT  over  the  Bay  of  Biscay  on 
August  20  1978  showed  a  number  of  internal  features  which 

were  generally  in  good  agreement  with  analysed  synoptic  and 
climatological  data  and  could  have  alternatively  be®n  used 
to  infer  valuable  sea-truth  in  operational,  real-time 
sampling. 

Wave  packets  propagating  offshore  from  the  shelfbreak 
had  an  inter-packet  spacing  on  tie  order  of  23KM  corre¬ 
sponding  *0  a  semi-diurnal  generation  at  the  main  canyons 
indenting  the  shelfbreak  at  low  oarotropic  tide  (maximum 
ebbing  surface  flow)  .  Speeds  computed  from  observed  density 
profiles  were  0.55M/S  ♦/-  0.5.  over  the  slope,  which 

compares  well  with  the  wave  speed  of  3.51M/S  calculated  for 
those  packets  from  the  SAR  image.  These  high  frequency 
internal  waves  have  a  period  of  the  order  of  70MIN. 
Assuming  these  waves  became  linear  far  enough  from  their 
generation  points,  a  mean  depth  of  the  mixed  lay«=r  of  20M 
was  inferred  and  spatial  variations  ascertained  from  clima¬ 
tological  data  and  observed  group  spacing. 

Drastic  spatial  variations  in  tae  wavelength  observed  in 
a  given  group  and  between  different  ones  were  analysed,  and 
assumed  to  be  due  to  current  shear  and  horizontal  density 
gradients.  They  coincided  with  a  sharp  thermal  front  at  the 
shelfbreak  and  an  alongshore  surface  gaostrophic  jet.  A 
line  of  surface  enhancement  of  Bragg  waves  corresponded  to 
this  front,  which  was  analysed  by  quasi-synoptic  X3T 
sampling.  This  front  is  generally  present  in  summer  condi¬ 
tions  when  a  corresponding  upwelling  breaks  the  surface. 
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Nonlinear  descriptions  of  *-he  internal  waves  in  each 
internal  wave  packet  yielded  estimates  of  their  amplitudes, 
probably  as  large  as  15M  on  the  continental  slope. 

A  tentative  application  to  underwater  acoustics  shows, 
under  severe  simplifying  assumptions,  that  detection  ranges 
in  such  a  situation  are  essentially  diraction-dependen*  and 
that  sonar  search  patterns  (or,  conversely,  submarine  escape 
routes)  could  be  greatly  enhanced  if  these  internal  features 
could  be  sampLed  on  a  regular  basis  and  taken  into  account. 
Adequate  operational  XBT  sampling  could  give  some  in-situ 
information  on  this  variability.  It  is  recommended  that,  in 
areas  where  internal  waves  are  known  to  occur,  consecutive 
launches  of  XBT's  be  made  at  a  time  (or  distance)  interval, 
corresponding  to  half  a  period  (wavelength) ,  thus  giving  an 
approximation  of  the  internal  wave  amplitude  and  conse¬ 
quently  of  the  variations  in  detection  ranges  that  can  be 
expected.  Internal  wave  climatology  developed  using  SAH 
images  would  be  helpful  in  designing  a  sampling  plan. 

The  following  ocean  features  of  operational  interest  can 
be  inferred  from  the  SAR  image  studied  here: 

a)  An  area  of  intense  internal  activity  is  clearly  visible 
on  the  continental  slope  in  the  form  of  short  internal  wave 
groups  propagating  from  fcha  heads  of  the  aajoc  canyons  and 
generated  approximately  at  the  time  of  low  tide  at  Brest. 
These  perturbations  probably  correspond  with  a  larger  scale 
widening  (weakening)  of  the  thermocline  and  periodic  mixing 
processes.  The  possible  operational  consequences  are 
(1)  drastic  changes  in  detection  ranges  due  to  a  varying 
thermal  structure  in  time  and  space  and  a  probable  increase 
in  ambient  noise  level  due  to  bottom  stress  and  turbidity 
currents,  and  (2)  sudden  modifications  in  buoyancy  that  may 
affect  underwater  vehicles  (e.g.,  submarines,  torpedoes)  . 

b)  A  well-defined  front,  located  at  the  upper  oart  of  the 
sheifbreak,  is  apparent  on  the  SAR  image  as  a  strong 
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localised  enhancement  of  the  high  frequency  waves.  It  corre¬ 
sponds  to  the  well-known  shelfbreak  upwelliag  in  the  summer 
in  the  area,  which  has  been  demonstrated  by  synoptic  XBT 
sampling.  Iasonif ication  through  such  a  front  is  essen¬ 
tially  dependent  on  the  relative  positions  of  source  and 
receiver. 

c)  In  the  open  sea,  several  fields  of  internal  waves  are 
organised  in  groups  spreading  from  the  shelf.  The  general 
direction  of  the  crests  is  parallel  to  the  isobaths  and 
creates  anisotropies  in  sound  speed  structure  and  inscnifi- 
cation  volume. 

d)  Other  features,  not  considered  here,  have  eventual 
tactical  importance:  (1)the  localisation  of  areas  of  calms 
and  therefore  mirror  sea  surface  for  visual  search  in  ASH 
and  SAR  (Search  and  Rescue)  missions  and  site  selection  for 
weapon  launch,  and  (2)  the  coastal  refraction  of  the  swell 
giving  a  possible  surf  hindcasting  capability  for  mine¬ 
sweeping  and  amphibious  landing  operations. 

The  coming  SIR-B  missions  will  cover  the  Bay  of  Biscay 
and  consequently  provide  some  information  on  the  time  vari¬ 
ability  of  the  characteristics  appearing  in  this  single  SAR 
image.  Combined  international  effort  to  increase  synoptic 
in-situ  data  acquisition  at  this  time  should  be  eminently 
ben ef icial. 

Spaceborne  imaging  radars  will  supply  new  capabilities 
for  understanding  and  perhaps  even  continually  sampling  the 
evolution  of  the  internal  ocean.  These  capabilities  may 
dramatically  change  the  parameters  of  undersea  warfare. 
Every  effort  should  therefore  be  made  in  allied  navies  to 
understand  and  master  these  consequences,  from  both  scien¬ 
tific  and  operational  points  of  view,  to  make  sure  that  no 
dangerous  and  destabilizing  consequences  could  result  from 
their  universal  exploitation. 
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